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Abetract 

On  the  basis  of  the  equation  of  heat  conduction,  a  mathematical  formula 
is  obtained  which  gives  the  time  constant  of  the  pneumatic  cell  as  a  function 
of  the  depth  of  the  cell.  The  numerical  calculation  of  the  time  constant  is 
carried  out  with  this  formula  assuming  plausible  values  of  various  physical 
constants. 

From  these  results  the  temperature  rise  in  the  cell  is  calculated  under  the 
assumption  that  the  cell  has  only  one  time  constant.  Thus  it  is  ascertained 
that  there  is  an  optimum  depth  of  the  cell  which  gives  the  maximum  temper¬ 
ature  rise  for  a  given  choiring  frequency. 

The  displacement  of  a  flexible  mirror,  which  is  a  measure  of  the  sensi¬ 
tivity,  is  also  estimated  using  the  equation  of  state  of  the  gas  and  the  tem¬ 
perature  rise  in  the  cell  already  obtained,  and  it  is  found  that,  under  ordinary 
conditions,  the  sensitivity  shows  a  maximum  at  a  certain  value  of  the  depth 
of  the  cell,  which  is  fairly  larger  than  the  value  corresponding  to  the  maximum 
temperature  rise.  % 

The  pressure  deptndence-  of  the  sensitivity  is  derived  from  the  combination 
of  the  various  formuj^  obtained,  and  it  is  shown  that  the  cell  has  its  maximum 
sensitivity  at  a  pre^ure  between  0.1  and  1  atm. 

Lastly,  the  signal-to-noise  ratio  of  the  cell  is  discussed  assuming  that  there 
is  no  noise  other  than  that  caused  by  temperature  fluctuation,  and  the  minimum 
detectable  power  is  also  estimated. 


1.  Introduction 

The  pneumatic  cell  invented  by  M.  J,  E.  Golay",  being  an  excellent  infra-red 
detector,  is  now  often  used  in  various  fields  of  infra-red  work,  especially  in  the 
study  of  the  far  infra-red.  Unfortunately,  very  few  data  on  its  performance  are 
available  and  one  find  difficulties  in  constructing  the  cell,  although  Golay*’  has 
established  the  fundamental  principle  of  construction. 

In  this  paper,  the  time  constant  and  the  sensitivity  of  the  cell  are  calculated 
by  directly  solving  the  equation  of  heat  conduction  to  discuss  the  important  char¬ 
acteristics  of  the  cell  and  to  give  the  detailed  data  on  its  performance. 

1)  H.  H.  Zahl  and  M.  J.  E.  Golay:  Rev.  Sci.  Instr.  17  (1946)  551. 

M.  J.  E.  Goay;  ibid,  18  (1947)  347;  18  (1947)  357. 

2)  M.  J.  E.  Golay;  Rev.  Sd.  Instr.,  20  (1949)  811. 
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2.  A  Rough  Estimation  of  the  Time  Constant 

We  shall  first  calculate  the  time  constant  of  the  cell  without  the  heat  absorbing 
film.  (Fig.  1) 

As  will  be  seen  later,  the  result  thus 
obtained  proves  to  be  a  fairly  good  approxi¬ 
mation.  For  simplicity  the  effect  of  the 
narrow  duct  leading  to  the  flexible  mirror 
as  well  as  of  the  pneumatic  leak  is  neglected 
for  the  calculation. 

Consider  a  cylindrical  pneumatic  cell, 
and  let  its  center  be  the  origin  of  co- 
odinates,  its  axis  the  z-axis.  (we  shall  use 
hereafter  the  cyrindrical  co-odinates) 

The  equation  for  the  heat  conduction 
*^‘8-  1  of  the  gas  is  given  by 

.=kJ^(jT„) 

where  JTg  is  the  temperature  rise  of  the  gas  above  the  ambient,  k  is  KgIpgCg  and 
Kg,  pg  and  Cg  are  the  thermal  conductivity,  the  density  and  the  heat  capacity  per 
unit  mass  of  the  gas  respectively. 

The  boundary  conditions  in  this  case  are 

JT,=0  when  z=±d 

JTg=0  when  r=a 

where  a  and  d  are  the  radius  and  the  half  of  the  depth  of  the  cylinder  respec¬ 
tively.  As  is  well  known,  the  solution  of  the  equation  is  given  by 

i 

sin  ^  (z-Hd)(A„. «, ,  cos  nO+B^, »,  n  sin  nO) 

where  J»  is  the  Bessel  function  of  the  n-th  order,  p  is  determined  by  /«(/<•  a) =0, 
and  A’s  and  B’s  are  the  constants  determined  by  the  initial  condition. 

This  solution  shows  that  the  pneumatic  cell  has  an  infinite  number  of  time 
constants  (namely  the  reciprocals  of  mV/4<f))  unlike  other  infra-red  detectors 
such  as  thermocouples  or  bolometers. 

In  practice,  however,  when  the  chopped  beam  method  followed  by  the  elec¬ 
tronic  amplification  is  utilized,  the  largest  time  constant  is  the  most  important 
so  far  as  the  chopping  frequency  is  sufficiently  low  and  the  other  time  constants 
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are  much  smaller  than  the  largest.  As  this  condition  is  nearly  satisfied  in  the 
present  case,  we  may  regard  the  pneumatic  cell  as  a  detector  with  only  one  time 
constant  which  corresponds  to  the  lai^est  one. 

Thus  the  time  constant  r  is  given  by 

r=l/*(A**oi+T*/4<?),  fM>ia=2.m  (2) 

For  example,  if  r=10"*sec.  and  a=0. 15cm,  d  is  equal  to  9. 8x10"* cm  for 
the  air-filled  cell.  Hence  we  see  that,  for  the  pneumatic  cell  used  with  the  chopped 
beam  method,  the  depth  of  the  cell  must  be  of  the  order  of  1  mm. 

For  the  “flat  cell”  for  which  a>d,  the  equation  (2)  becomes 


3.  Calculation  of  the  Time  Constant 

We  now  calculate  the  time  constant  of  the  cell  with  a  heat  absorbing  film 
within  it.  Let  this  film  coincide  with  the  ay-plane  in  Fig.  1. 

According  to  the  usual  method  of  solution  we  assume  dT,=U{t)  R(r)  6(0)  Z(z), 
where  U,  R,  6  and  Z  are  the  functions  of  t,  r,  0  and  z  respectively.  From  physical 
considerations  JT,  must  be  axially  symmetric,  and  consequently  tP6ld0^=O.  Thus 
we  have 


1 

dU 

=ki-^ 

(PR 

+  L 

1  dR  . 

1 

(PZ  \ 

u 

dt 

\r 

dr* 

^  R 

r  dr 

dz*  ) 

Put  this  expression  equal  to  —k(a*+^),  and  let  (\IZ)(cPZIdz^)=—^,  where  a 
and  p  are  constants.  Then  we  have 

1  ^R  ,  1  dR  __  , 

R  dr  ^  R  r  dr  " 

1  «  .  . 

-z  </«■ 

The  particular  solution  of  the  problem  is  given  by 

JT„=exp{— *(a*-fy3*)/}y»(ar)(ci  cos  pz-\-ci  sin  pz)  ( 3 ) 

where  Ci  and  Cz  are  constants. 

Boundary  conditions  in  this  case  are 


o 

II 

when 

z=d 

(4). 

JT,=0  when  r=a 

(5) 

u 

II 

u 

when 

z=0 

(6). 

II 

d(dT,)  \ 
dz  Jm-o 

(7) 
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where  AT f  is  the  temperature  rise  of  the  absorbing  film,  c/  is  the  heat  capacity 
of  the  film  per  unit  area  and  kt  is  the  radiation  cooling  constant  per  unit  area. 
(Note  that  cg  and  kt  have  different  dimensions) 

At  first  sight  it  seems  strange  that  the  condition  JTg=0  when  z=—d  is  omitted, 
but  addition  of  this  condition  gives  rise  to  a  contradiction.  This  is  because  the 
absorbing  film  has  in  reality  a  finite  thickness  and,  strictly  speaking,  AT,  is  not 
zero  when  z=—d  if  the  origin  of  the  co-ordinates  is  on  the  right  side  of  the  film. 
To  avoid  this  difficulty,  we  consider  only  the  right  half  of  the  cell  and  double 
both  the  radiation  and  conduction  losses  as  shown  in  the  condition  (7). 

From  the  condition  (4)  and  (5),  AT,  given  by  (3)  becomes 

JTa=const.  exp[—k{/i,^+ sin ^d—z)  ( 8 ) 

where  //,  is  determined  by  /o(/i.a)=0. 

The  condition  (6)  gives 

J  7/ = const.  sin(,5</)exp  { -  -I-  /3*)f)/o(/<tr)  ( 9 ) 

Substituting  the  expression  (8)  in  (7),  we  have 

=  -2crJr/-2*,  const.  ^•exp{-*(/«.*-f ^*)f}/o(//.r)  (10) 

The  equations  (9)  and  (10)  give 

=  -  {2irr+ cot(^<0}  AT, 

Therefore 


AT, =const.'  exp{-^^i±2f[£|E2iM_/} 

This  equation  must  be  identical  with  (9).  Consequently  we  have 
2fr + 2xg^  cot(,3rf) = -I-  j3*)c/ 

This  equation  determines  the  constant  ^9.  Since  the  cotangent  is  a  periodic 
function,  there  are  an  infinite  number  of  ;3’s.  The  general  solution  for  AT,  is 
given  by  summing  up  the  particular  solutions  such  as  (3)  with  various  values  of 
a  and 

Writing  r=  {^(a/.*-I-^*)}“‘,  we  have 

arr+2r,  V  (l/*r)-/*.*COt(rfV  (l/kr)- fi,*  )=c,It 
or  solving  for  d,  we  have 


d=  —  ^  .  —  arccot ( 

V  (l/jtr)-/ii*  V  ar,V(l/ih-)-«i*  / 


(11) 


From  the  same  reason  described  in  the  preceeding  section,  we  are  concerned 
with  the  largest  time  constant  only.  Hence  we  may  put  and  take  the 

principal  value  of  the  arccotangent. 
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For  the  flat  cell  for  which  a>rf,  /ii=0  and  the  equation  (11)  becomes 

d=Vjbarccot(-^£^/^^v^)  (12) 

These  equations  can  be  used  to  calculate  d  from  r,  and  allows  us  to  design 
the  cell  suitable  for  our  purpose. 

We  note  here  that  (1/^r)— ;/i*  must  not  be  negative  in  (11),  and  thus  there 
exists  no  cell  of  which  the  time  constant  exdeeds  Vftx^k. 


4.  Approximate  Formulas  for  the  Time  Constant 


In  the  case  of  the  flat  cell  we  can  derive  relatively  simple  formulas  for  the 
time  constant  from  the  equation  (12). 

Case  1.  t<C/I2kt 

Expanding  the  cotangent,  we  have 


C//t— Zr, 


Hence 


c/ 


2«r  +  (2<,/d) 


(13) 


The  time  constant  is  the  heat  capacity  of  the  film  divided  by  the  total  cooling 
constant.  The  heat  capacity  of  the  gas  is  in  this  case  so  small  that  it  may  be 
neglected. 

Case  2.  T=Cfl2KT 

Expanding  the  arccotangent  in  the  neighbourhood  of  zero  of  the  argument, 
we  have 


Solving  for  r,  we  obtain 

c/-f(2g^/^) 

Because  r  is  approximately  equal  to  ffVjkr/Z  in  this  case,  we  have 


_ c/-k-(2K,dlk) 

2kt+ (xl2y(2K,ld) 


(14) 


If  we  note  here  that  k  is  kqIp^q,  2r«<f/A  is  equal  to  2dp^g,  which  is  the  heat 
capacity  of  the  gas  per  unit  cross  section.  On  the  other  hand,  (xl2y(2K,ld)  re¬ 
presents  the  “effective”  cooling  caused  by  the  conduction  of  heat  through  the 
gas.  Thus,  the  physical  meaning  of  the  equation  (14)  is  that  the  time  constant 
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of  the  cell  is  the  total  heat  capacity  (the  heat  capacity  of  the  film  and  the  gas) 
divided  by  the  total  cooling  constant  (the  radiation  cooling  plus  the  conduction 
cooling). 

Furthermore,  it  is  clear  that  the  equation  (14)  is  reduced  to  the  equation  (2)’ 
in  the  section  2,  if  CfK2Kfdlk  and  2/cr<(ff/2)*(2/r,/</). 

Case  3.  r>  C//2kt 

Using  the  formula  arccot  - )  ^<0,  we  have 

Hence 

-  %  («) 

This  equation  shows  that  the  time  constant  is  chiefly  determined  by  the  heat 
capacity  of  the  gas  and  the  conduction  through  the  gas,  but  the  radiation  cooling 
is  to  be  taken  into  account  as  a  correction  term. 

As  these  formulas  are  very  simple  and  have  clear  physical  meaning,  they  are 
often  useful  to  anticipate  the  main  characteristics  of  the  pneumatic  cell. 

5.  Numerical  Calculation  of  the  Time  Constant 

In  order  to  evaluate  accurately  the  time  constant  of  the  pneumatic  cell,  we 
must  use  the  formula  (11)  or  (12).  For  this  purpose,  it  is  necessary  to  know  the 
physical  constants  which  appear  in  these  formulas. 

Since  no  definite  values  concerning  the  physical  constants  of  collodion  are 
available,  we  assume  that  the  specific  heat  is  0.4cal/g  deg.,  the  density  is  1. 5g/cm*, 
and  furthermore  the  thickness  is  400  A.  Then,  the  the  capacity  of  collodion  film 
per  unit  area  is  2. 4x  10"®  cal/deg-cm*.  In  the  present  experiment  we  used  antimony 
as  heat  absorbing  material.  The  data  for  antimony  are  easily  found  in  any 
physical  table.  If  we  assume  the  thickness  of  antimony  film  is  200  A,  the  heat 
capacity  of  antimony  becomes  6. 65xl0"’'cal/deg-cm^. 

Thus,  for  the  total  heat  capacity  of  the  absorbing  film  per  unit  area  we  have 
c/=3.  Oe  X  10“®  cal/deg-cm* 

As  for  the  gas  used  in  the  cell,  xenon  is  the  most  ideal  because  it  has  the 
smallest  thermal  conductivity  among  all  kinds  of  gases  suitable  for  the  pneumatic 
cell.  It  is,  however,  interesting  to  calculate  the  data  for  the  air-filled  cell  and  to 
compare  them  with  those  of  xenon. 

Since  the  gas  has  two  kinds  of  heat  capacity,  namely,  the  heat  capacity  at 
constant  volume  and  that  at  constant  pressure,  there  is  some  ambiguity  as  to 
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which  heat  capacity  is  adequate  for  our  purpose.  In  the  present  calculation  the 
capacity  at  constant  pressure  is  employed,  because  in  the  actual  cell  the  gas 
expands  to  some  extent  when  it  is  warmed  by  the  incident  radiation.  There  is, 
however,  little  difference  if  the  capacity  at  constant  volume  is  assumed. 

Thus,  for  xenon  we  have 

Ixl0“®cal/cm  sec.  deg. 

— =4.96  10“*cm*/sec.  (at  1  atm.) 

PfCg 

and  for  air  we  have 

<r»=5. 65x10"®  cal/cm  sec.  deg. 
kg=~^ — =1. 95x10"' cmVsec.  (at  1  atm.) 

pgCg 

The  radiation  colling  constant  can  be  easily  evaluated  using  the  Stefan-Boltz- 
mann’s  law  and  we  have 

cr=7. 35x10"®  cal/cm*  sec.  deg. 

where  we  have  assumed  that  the  emissivity  of  the  absorbing  film  is  0.5. 

Once  these  constants  are  determined, 
the  dependence  of  the  time  constant 
on  the  distance  between  the  absorbing 
film  and  the  wall  (d)  can  be  easily  cal¬ 
culated  numerically  from  (11)  and  (12). 

The  results  are  shown  in  Fig.  2. 

Experimental  results  so  far  obtained 
seem  to  support  this  calculation. 

It  is  seen  from  this  figure  that, 
under  similar  conditions,  the  xenon- 
filled  cell  has  a  larger  time  constant 
than  the  air-filled  cell.  But  this  does 
not  necessarily  mean  the  disadvantage 
of  the  former.  On  the  contrary,  on 
account  of  low  thermal  conductivity, 
the  temperature  rise  in  the  xenon  cell 
becomes  much  higher,  and  as  will  be 
shown  in  the  following  sections,  the  xenon  cell  is  superior  to  the  air-filled  one. 

6.  The  Mean  Value  of  the  Temperature  Rise  in  the  Cell 

We  now  proceed  to  calculate  the  mean  value  of  the  temperature  rise  in  the 
cell  which  is  of  primary  importance  in  the  problem  of  sensitivity,  although  the 


Fig.  2 
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sensitivity  depends  not  only  on  the  temperature  rise  but  also  on  other  factors  as 
will  be  seen  in  the  next  section. 

First  we  calculate  the  temperature  rise  for  constant  incident  radiation.  The 
equation  to  be  solved  is 

with  the  boundary  conditions 

JTf=0  when  z=d,  J7’„=0  when  r=a, 

JTt=JT/  when  2=0 

Q=2KrdTf-2K,(  ) 

\  OZ  /f-0 

where  Q  is  the  power  actually  absorbed  per  unit  area.  (Q  is  in  general  a  function 
of  r) 

Using  a  similar  procedure  as  in  the  section  3,  we  have 
JTt=A’Jo{ttr)  sinh  {ft(d—z)} 
where  ft  is  given  by  /o(/<«)=0,  and  ^4  is  a  constant. 

Also  we  have 

JT /=AJo{fir)  sinh  (fid) 

Then,  the  last  of  the  boundary  conditions  becomes 

^ _  Q = TxrAJfitir)  sinh  {ptd) + 2KgAJ(iiir)fi  cosh  (/«/) 

For  simplicity  we  assume  that  the  distribution  of  the  incident  radiation  is  not 
uniform  but  is  given  by  Q=QoJ<^fir).  (Otherwise  Q  must  be  expanded  in  a  series 
of  Bessel  functions) 

Thus  we  obtain 

_ Qo _ 

2kt  sinh  {/id) + 2k, ft  •  cosh  (/«/) 

and  '  t 

The  mean  value  of  the  temperature  rise  in  the  cell  is  given  by 

+ 2«,  •  //  •  cosh  {(id) 

X  jyo(/tr)-2rr- rfrj^sinh{/i(</-2)}</2 

=  Oo  co^  C/ad)-!  _  ^  _  (16) 

[la  (id  sinh  ijid)  2«r+2«,/:i-coth  {(id) 
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For  the  flat  cell  we  have 


In  these  formulas,  the  suffix  0  of  JT,  means  that  the  chopping  frequency  is  zero. 
In  the  case  of  the  flat  cell,  if  the  incident  radiation  is  uniform  over  the  ab¬ 
sorbing  film,  we  have  by  a  simple  calculation 


1 

2xr  +  2xf  Id 


This  is  different  from  (17)  by  a  factor  2Mua)l(ui.  The  difference  arises  from 
the  fact  that  the  incident  radiation  is  assumed  non-uniform  in  the  above  calculation. 


If  we  assume  that  the  pneumatic  cell  has  only  one  time  constant  given  in  the 
preceeding  section,  the  estimation  of  response  for  chopped  incident  radiation  can 
be  carried  out  easily  as  in  the  case  of  thermocouples  or  bolometers.  The  result  is 


(JT,).=(jr,)o 


where  v  is  the  chopping  frequency. 

(JT»)o  is  an  increasing  function  of  d  for  small  d,  whereas  l/(l-l-(2rvT)*)  is  a 
decreasing  function  because  r  increases 
with  d.  Consequently  ( JT »),  becomes 
maximum  at  a  certain  yalue  of  d.  Thus  po 
we  see  there  is  an  optimum  depth  of  ^ 
the  cell  for  a  given  chopping  frequency. 

In  Fig.  3  (JT ,).  given  by  (18)  is  ^ 
plotted  against  d  assuming  Qo=10~^  H 

cal/cm*  and  v=10c/s.  The  thickness 
of  the  absorbing  film,  the  physical  con- 
stants,  etc.  are  the  same  as  used  in 
the  section  5. 

This  figure  shows  that  the  xenon 
cells  are  decidedly  more  advantageous  — '  ^  (cm) 

than  the  air  ones  so  far  as  the  mean  ^ 


iiiiiifim 
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temperature  rise  is  concerned.  Furthermore,  the  loss  due  to  finite  diameter  of 
the  cell  is  not  so  serious  as  was  first  expected,  especially  for  the  xenon  cells.  Thus 
we  see  that  the  increase  of  diameter  makes  little  gain  in  the  mean  temperature  rise 
as  far  as  the  distribution  of  the  incident  radiation  over  the  absorbing  film  is  similar. 

Unfortunately,  these  calculations  cannot  be  verified  experimentally,  because  it 
is  difficult  to  construct  similar  cells  except  the  depth,  and  furthermore  there  is  no 
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means  to  measure  directly  this  temperature  rise  of  the  gas. 

We  must  note  here  that  the  equation  (18)  is  not  valid  for  all  values  of  d.  The 
wave-length  of  the  heat  wave  generated  by  periodical  irradiation  becomes  shorter 
with  increasing  frequency,  and  when  it  finally  becomes  smaller  than  the  depth  of 
the  cell,  maxima  and  minima  of  temperature  appear  in  the  cell  and  it  becomes 
impossible  to  regard  the  cell  as  a  detector  with  only  one  time  constant.  In  the 
calculation  above,  where  the  chopping  frequency  is  10  c/s,  the  wave-length  of  the 
heat  wave  is  about  5  mm  for  air  and  2  mm  for  xenon.  Thus  we  see  that  the 
equation  (18)  is  a  good  approximation  at  least  in  the  neighbourhood  of  the  maximum 
in  Fig.  3,  though  near  the  right  ends  the  curves  do  not  correctly  represent  the 
real  behavious  of  the  cell. 


7.  Sensitivity  of  the  Cell 

The  pneumatic  cell  is  the  more  sensitive  the  larger  the  displacement  of  the 
flexible  mirror  is,  provided  that  the  intensity  of  the  incident  radiation  is  the  same. 
Hence  we  may  take  this  displacement  as  the  measure  of  sensitivity. 

We  now  consider  how  this  displacement  depends  on  various  parameters  of 
the  cell. 

The  equation  of  state  of  the  gas  is 

pV=NRTg 

where  p,  V  and  T,  are  the  pressure,  the  volume  and  the  temperature  of  the  gas 
respectively.  N  is  the  number  of  moles  and  R  is  the  gas  constant. 

If  the  mean  temperature  rise  of  the  gas  is  Jl'g,  we  have 

J  V-\-  V^=NRWg  =  ^Wg 

Let  the  area  of  the  flexible  mirror  be  A  and  its  effective  displacement  be  Jx, 
then  we  have  JV=A-Jx. 

Let  the  constant  of  the  restoring  force  of  the  mirror  be  a,  and  assume  that 
the  space  outside  the  mirror  is  so  large  that  its  pressure  does  not  change  ap¬ 
preciably.  Then  we  have  A-Jp=a-Jx. 

From  these  relations  it  follows  that 

Jx—  ^P'd'f  t  QQX 

na+pA*lV) 

Strictly,  the  above  discussion  is  valid  only  for  constant  incident  radiation,  but 
the  natural  frequency  of  the  mirror  is  so  high  when  compiared  with  the  usual 
chopping  frequency  that  these  formulas  are  practically  good  approximations.  (The 
natural  frequency  can  be  evaluated  using  the  data  for  surface  tension  obtained 
by  (jolay.) 


Properties  of  Pneumatic  Infra-red  Detector  (I) 


77 


From  the  expression  (19)  we  can  calculate  the  sensitivity  of  the  cell  as  follows. 
Case  1.  Flat  cell 

In  this  case  the  diameter  and  consequently  the  volume  of  the  cell  are  so  large 
that  PA*IV  is  neglibible  as  compared  with  a.  Then  we  have 

7a 


If  A,  p,  T  and  a  are  fixed,  Jx  is  proportional  to  JT,. 

Case  2.  Cell  with  finite  diameter 

In  this  case  the  term  pA*lV  becomes  important  as  is  seen  by  the  following 
estimations. 

The  restoring  force  consists  of  the  elastic  force  and  the  surface  tension  of 
the  him.  As  the  flexible  mirror  is  made  of  extremely  thin  membrane,  the  elastic 
force  may  be  negligible  when  compared  with  the  surface  tension.  If  this  is  the 
case,  the  constant  of  the  restoring  force  a  proves  to  be  simply  AkP  by  an  ele¬ 
mentary  calculation,  where  f  is  the  constant  of  the  surface  tension.  According 
to  Golay  T  is  of  the  order  of  10*  dynes/cm,  whence  AkP  is  of  the  order  of  10* 
dynes/cm.  On  the  other  hand,  PA^jV  is  about  2. 4x10*  dynes/cm  at  1  atm.,  assum¬ 
ing  that  a=rf=0. 15cm  and  A=7xl0“*cm*.  Thus  we  see  that  AnP  is  generally 
fairly  smaller  than  PA*IV. 

Then  we  have  for  the  displacement  of  the  mirror 


Jx= 


V-dT, 

TA 


We  find  that  in  this  case  the  sensitivity  depends  upon  the  volume  of  the 
cell,  and  in  particular  it  depends  upon  the  depth  of  the  cell  if  the  diameter  of  the 
cell  is  fixed.  On  account  of  this  dependence,  the  curve  of  Jx,  plotted  against  d, 
becomes  very  different  from  that  of  JT ,,  the  maximum  being  displaced  tow^ds 
a  larger  value  of  d. 

In  Fig.  4  are  shown  the  sensitivity  ciuwes  calculated  from  the  above  expres¬ 
sions.  For  the  flat  cells  the  positions  of  maxima  agree  with  those  of  JT ,,  whereas 
for  the  cells  with  finite  diameters  the  maxima  occur  at  about  1. 5  mm,  though  they 
are  not  so  distinct. 

We  must  remember  that,  as  was  pointed  out  in  the  preceeding  section,  the 
curves  in  Fig.  3  and  accordingly  the  curves  in  Fig.  4  do  not  represent  the  real 
behaviour  of  the  cell  near  the  right  ends.  Therefore,  it  is  a  question  whether  the 
maxima  of  the  curves  in  Fig.  4  correspond  precisely  to  the  real  maxixa  of  the 
sensitivity.  Perhaps  the  sensitivity  will  fall  more  rapidly  than  these  curves  when 
a  exceeds  1  or  2  mm,  because  the  temperature  distribution  begins  to  show  a 


78 


Kunio  Yosihara 


minimum  near  the  wall  of  the  cell,  and 
as  a  result  of  this,  the  positions  of  these 
mixima  will  be  somewhat  displaced  to 
the  left.  But  a  more  exact  calculation  is 
necessary  to  get  correct  sensitivity  curves. 

Usually  we  employ  the  pneumatic 
cell  instead  of  thermocouples  when  the 
incident  beam  has  so  large  a  cross  sec¬ 
tion  that  the  thermocouple  cannot  utilize 
all  of  the  incident  power  on  account  of 
its  small  receiving  area,  as  is  the  case 
with  the  far  infra-red  spectrometer.  In 
such  a  case  it  seems  advantageous  to 
make  the  area  of  the  absorbing  him 
just  equal  to  the  cross  section  of  the 
incident  beam,  because  the  sensitivity 
increases  remarkably  with  increase  of 
the  diameter  of  the  cell  as  shown  in 
Fig.  4.  We  must  be  careful,  however, 
to  distinguish  the  sensitivity  from  the 
signal-to-noise  ratio.  It  is  because  the  factor  ApjT  {a+pA^jV)  in  (19)  becomes 
large  that  the  flat  cell  has  a  greater  sensitivity  than  that  of  the  cell  with  flnite 
diameter.  But,  for  the  same  reason,  the  fluctuation  of  temperature  in  the  cell 
(temperature  noise  etc.)  gives  rise  to  a  large  random  displacement  of  the  flexible 
mirror.  Hence,  increase  of  the  factor  ApiT{a-\-pA*IV)  does  not  mean  the  improve¬ 
ment  of  the  signal-to-noise  ratio.  This  ratio  is  limited  in  principle  by  the  tem¬ 
perature  fluctuation  in  the  cell  itself,  if  there  is  no  other  kind  of  fluctuation. 

The  signal-to-noise  ratio  of  the  cell  will  be  discussed  later. 

> 

8.  Effect  of  Gas  Pressure  on  the  Sensitivity 

It  is  interesting  to  investigate  the  dependence  of  the  sensitivity  on  the  gas 
pressure  in  the  cell.  Indeed,  this  dependence  is  generally  very  much  complicated 
and  it  is  almost  impossble  to  estimate  it  exactly,  but  some  qualitative  conclusions 
can  be  obtained  by  applying  the  formulas  already  derived. 

Let  the  chopping  frequency  v  and  the  depth  of  the  cell  be  fixed,  and  let  the 
gas  pressure  only  be  variable.  For  simplicity  we  shall  use  the  expression  (14)  for 
the  time  constant  and  the  expression  (17)  for  (JT ,)%,  even  when  the  diameter  of 
the  cell  is  not’ very  large. 
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We  must  remember  that  the  heat  conductivity  of  the  gas  does  not  depend 
on  the  pressure  so  far  as  the  pressure  is  higher  than  the  range  in  which  molecular 
flow  takes  place. 

( 1 )  Cell  with  finite  diameter 

Low  frequency  operation 

When  the  chopping  frequency  is  sufficiently  low,  namely  (2.Tvr)*<l,  we  have 
from  (18)  and  (19) 

Zcr+W/d  T{a-^pA*lV) 

The  time  constant  does  not  appear  in  this  expression.  When  the  gas  pressure 
P  is  higher  than  1  atm.,  a  is  small  as  compared  with  pA^jV  and  we  have 


2«,+2c,/d  TA 

As  is  independent  of  p  in  this  range,  Jx  must  be  nearly  constant. 

Since  a  is  about  one  tenth  of  pA^j  V,  a  cannot  be  neglected  when  p  approaches 
10"‘atm.,  and  then  Jx  begins  to  decrease. 

When  the  pressure  is  considerably  below  10"‘atm.  but  is  greater  than  the 
range  of  molecular  flow,  the  expression  (20)  becomes 

Jr-  1 

2er-i-2K,ld^ 

Thus  we  see  that  Jx  decreases  linearly  with  the  decreasing  pressure  in  this 


region. 

As  the  diameter  of  the  cell  is  of 
the  order  of  one  millimeter,  the  mol¬ 
ecular  flow  begins  to  occur  at  about 
10"*  atm.,  and  the  heat  conductivity 
K,  also  begins  to  decrease.  As  a  result 
of  this,  the  rate  of  decrease  of  Jx 
becomes  somewhat  smaller. 

If  the  pressure  is  reduced  further, 
the  heat  conductivity  becomes  so  small 
that  it  can  be  neglected  as  compared 
with  radiation  cooling.  Thus  we  have 


Ap 

Ta 


/O'*  id*  10*  10*  w*  to'  I.V  to 

- ►  Preaeure  (Atm.) 

Fig.  5 


Hence,  Jx  decreases  again  linearly  with  decreasing  pressure. 

The  curve  A  in  Fig.  5  indicates  such  a  variation  of  sensitivity  with  pressure 
as  a  whole. 
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High  frequency  operation 

In  this  case  (2ffvr)*>l,  and  it  follows  from  (18)  and  (19)  that 


ijr-  1 _ AP  „  1 

2rr+2r,ld  T(a+pA*IV)  t 

where  r  is  given  by 

_  Cf-\-  ZdCgpg 

at.+(r/2)*(at,/rf) 


according  to  (14).  In  this  expression,  2dpgC,  proves  to  be  about  ten  times  greater 
than  c/  at  latm.  if  d  is  about  1.5  mm.  Hence,  when  p  is  greater  than  10"*  atm., 
we  have  On  the  other  hand  we  have  aKpA^IV  as  in  the  previous  case. 

Thus  it  follows  that 


2rr+2r, Id  TA  p 

Since  tg  is  independent  of  p  in  this  range,  Jx  is  inversely  iHX>portional  to  the 
gas  pressure. 

When  10"*</><10"*  atm.,  a  and  c/  are  far  greater  than  pA^IV  and  2dCgpg 
respectively.  Then  we  have 


Jxoc 


1 

2Kr  +  2Kgld 


AP^l 

Ta  T 


T  = 


Cf _ 

2jCT  +  {i:l2Y&Kgld) 


In  this  range  of  gas  pressure.  Kg  is  independent  of  pressure  and  consequently 
Jx  is  proportional  to  p. 

At  a  pressure  below  10"‘atm.,  begins  to  decrease,  whence  \l{2Kr-\-2xgld)  as 
well  as  T  begins  to  increase.  Thus  Jx  decreases  more  slowly  as  the  pressure  is 
reduced.  At  much  lower  pressure,  becomes  so  small  that  we  have 


Jxoc 


1 

2Kr 


Ap  1 
la  r 


Hence  Jx  becomes  again  proportional  to  the  gas  pressure.  The  curve  B  shows 
the  variation  of  sensitivity  in  this  case. 

Intermediate  frequency  operation 

In  actual  cases,  (2;rj/r)*  has  a  value  comparable  with  1  and  the  pressure  de¬ 
pendence  of  sensitivity  becomes  complicated.  But  the  curve  representing  this 
dependence  must  lie  between  the  curves  A  and  B  in  Fig.  5.  Thus  the  broken 
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curve  such  as  C  may  be  used  as  a  rough  approximation.  We  may  expect  from 
the  general  tendency  of  this  curve  that  the  pneumatic  cell  has  a  maximum  sensi¬ 
tivity  at  a  pressure  between  10~*  and  1  atm.,  though  this  maximum  is  presumably 
not  very  conspicuous. 

(2)  Flat  cell 

In  the  above  discussion  we  have  assumed  that  the  diameter  of  the  cell  is  not 
very  large  but  large  enough  to  permit  the  application  of  the  expressions  (14)  and 
(17).  For  the  flat  cell  the  situation  is  fairly  different.  In  this  case  the  volume  V 
becomes  very  large  owing  to  the  large  diameter  even  if  the  thickness  of  the  cell  is 
small.  Hence  a  is  large  when  compared  with  pA*IV  in  (19)  and  we  have 

Jx=  (21) 

1  a 

Low  frequency  operation 

As  (2ffvT)*  is  by  far  smaller  than  1,  we  have  from  (21) 

Jxoc— JL__  AP. 

2cT-\-2Kfld  Ta 

We  can  see  from  this  expression  that  the  sensitivity  is  proportional  to  the 
pressure  near  latm.,  and  becomes  constant  when  the  pressure  is  so  large  that 
pA*lV  exceeds  a.  But  in  practice  such  a  high  pressure  cannot  be  attained. 

High  frequency  operation 

In  this  case  we  have 


Jxoo 


1 

2cr-|-2c«/</ 


AP  1 

Ta  T 


_ Cf-^2dc,p, 

‘  2«r-f(ff^)*(at,/d) 

From  the  same  reason  described  above,  C/  may  be  neglected  when  compared 
with  2dCgpg  if  the  pressure  is  comparable  with  1  atm.  or  greater.  Therefore* we 
have  as  before  roc^ocp.  As  a  result  of  this,  the  sensitivity  is  independent  of 
pressure  in  this  range.  It  is  not  until  the  pressure  is  lowered  below  10" ‘atm. 
that  the  sensitivity  begins  to  decrease. 

Intermediate  frequency  operation 

In  this  case  the  sensitivity  curve  can  be  anticipated  from  the  above  discussion. 
We  may  conclude  that  the  sensitivity  will  increase  with  increase  of  gas  pressure 
of  the  cell. 


9.  Signal  to-noise  Ratio  and  Minimam  Detectable  Power 


The  most  significant  criterion  for  the  infra-red  detector  is  the  signal-to-noise 
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ratio  rather  than  the  sensitivity  itself.  For  thermocouples  or  bolometers,  this  ratio 
is  chiefly  limited  by  the  Johnson  noise,  although  the  temperature  noise  becomes 
important  in  certain  conditions.  Unfortunately  the  origin  of  the  noise  other  than 
the  temperature  is  not  well  investigated  for  the  pneumatic  cell.  Indeed  Smith 
and  others’  pointed  out  that  the  damping  of  the  motion  of  the  flexible  mirror 
gives  rise  to  a  kind  of  undesirable  noise,  but  the  amount  of  this  noise  cannot 
be  estimated  numerically  because  the  nature  of  the  damping  are  not  clear. 
Hence  we  assume  for  the  present  that  there  is  no  noise  except  that  caused  by 
the  temperature  fluctuation,  and  discuss  the  results  derived  under  this  assumption. 

The  mean  temperature  fluctuation  JTnoiae  is  given  by” 


(JTaoiae)^  — 


(22) 


where  k  is  the  Boltzmann  constant,  k  is  the  total  cooling  constant,  C  is  the  total 
heat  capacity  and  Jv  is  the  frequency  bandwidth.  In  the  present  case  it  is  more 
reasonable  to  calculate  k  using  the  expression  k=CIt  than  to  calculate  it  directly, 
because  the  pneumatic  cell  is  assumed  to  have  only  one  time  constant. 

Thus  (22)  becomes 


(rlOJv 

l+(2ffvr)* 


Then  we  have  for  the  signal-to-noise  ratio  S/N 


S/iV=-  = 


1 

y/AknJv 


(23) 


Remembering  that  C=:ro’c/+2d(,-ro)’|0»c,,  S/N  can  be  evaluated  numerically 
with  this  expression.  The  results  are  shown  in  Fig.  6  in  which  S/N  is  plotted 
against  d,  assuming  a=0. 15cm,  Oo=10"’ cal/sec.  cm*  and  Jv=lc/s. 


It  is  seen  from  this  figure  that  there  is  an  optimum  depth  which  gives  the 
largest  signal-to-noise  ratio.  This  depth  is  somewhat  smaller  than  that  corre¬ 
sponding  to  the  maximum  sensitivity  but  is  larger  than  that  corresponding  to 
the  maximum  of  (JT »)„. 

In  the  case  of  the  cell  with  a  large  receiving  area,  however,  we  are  chiefly 
concerned  with  the  variation  of  S/N  with  the  area  of  the  absobing  film,  for  such 
a  cell  is  usually  employed  to  measure  a  beam  with  large  cross  section.  Let  the 
area  of  the  absorbing  film  be  E,  then  c  and  C  are  proportional  to  E,  and  accord- 


3)  R.  A.  Smith,  F.  E.  Jones  and  R.  P.  Chasmar;  The  Detection  and  Measurement  of 
Infra-red  Radiation.  Oxford  at  the  Qarendon  Press  (1957)  p.  278. 

4)  ibid.  p.  206. 
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ingly  the  mean  temperature  fluctuation 
is  inversely  proportional  to  On 

the  other  hand,  it  is  seen  from  Fig.  3 
that  ( JT j).  will  be  nearly  constant 
when  E  becomes  fairly  large.  Thus 
we  find  that  the  signal-to-noise  ratio 
is  approximately  proportional  to  the 
square  root  of  the  receiving  area  for 
the  flat  cell. 

The  effect  of  gas  pressure  on  the 
signal-to-noise  ratio  is  interesting.  The 
signal-to-noise  ratio  is  given  by  (23)  as 
before.  If  we  use  the  simplified  ex¬ 
pression  (14)  for  the  time  constant  r 
and  remember  that  the  numerator  of 
(14)  is  the  heat  capacity  per  unit  cross 
section,  we  obtain 

siN=  ,  - 

V4kt*Jv  fta 


{2KrH^I2y{2K^ld)Y<'E^i' 

2cr  +  2iCf/dr 


This  expression  does  not  contain  any  quantity  which  depends  upon  the  pressure 
so  far  as  the  pressure  is  higher  than  the  range  where  the  molecular  flow  occurs. 
Thus  the  change  of  pressure  has  no  effect  on  the  signal-to-noise  ratio  of  the  cell. 
Indeed,  use  of  more  exact  expression  of  r  may  show  some  depencence  of  SjN  on 
p,  but  we  may  conclude  that  this  depencence  is  very  small. 

ITie  minumum  detectable  power  can  be  calculated  easily  from  Fig.  6.  We 
assume  here  that  this  power  is  given  by  the  condition  S/iV=l.  If  we  use  the 
cell  in  the  optimum  condition  and  remember  that  the  absorptive  power  of  the 
receiving  film  is  0.5  and  the  distribution  of  the  incident  power  is  given  by  Qo 
x/o(//ir),  we  obtain  the  minimum  detectable  power  as  8. 6x10  “  watt  for  the  xenon 
cell  and  1. 6xl0~**  watt  for  the  air  cell,  provided  that  the  diameter  is  each  3  mm. 
The  value  for  xenon  cell  agrees  approximately  with  that  obtained  by  Golay. 
Therefore  the  minimum  detectable  power  of  the  pneumatic  cell  seems  to  be  limited 
chiefly  by  the  temperature  noise.  But  more  detailed  investigations  are  necessary 
about  this  problem. 


10.  Summary 


Starting  from  the  equation  of  heat  conduction,  we  are  able  to  calculate  the 
time  constant  and  the  sensitivity  of  the  pneumatic  cell,  and  using  these  results. 
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we  can  also  estimate  other  properties  such  as  the  signal-to-noise  ratio,  the  depen¬ 
dence  of  sensitivity  on  the  gas  pressure,  etc.,  so  far  as  the  cell  can  be  regarded 
as  a  detector  with  only  one  time  constant.  These  calculations  will  be  useful 
in  constructing  pneumatic  cells,  because  it  is  possible  to  anticipate  the  character¬ 
istic  properties  of  the  cells. 

The  numerical  results  depend  indeed  upon  the  assumption  about  the  thicknes 
of  the  collodion  film,  the  value  of  the  specific  heat  of  collodion,  the  value  of  the 
surface  tension  of  the  flexible  mirror,  etc.,  but  the  essential  part  of  the  discussion 
will  not  be  altered. 

Strictly  speaking,  the  pneumatic  cell  hzis  an  infinite  number  of  time  constants 
and  therefore  it  is  desirable  to  calculate  more  rigorously  the  response  to  the  alter¬ 
nating  incident  radiation.  Such  a  calculation  will  be  published  in  another  paper. 

Experimental  results  obtained  in  our  laboratory  will  be  also  published  in  the 
near  future,  but  we  would  add  here  that  recordings  of  many  infra-red  spectra 
were  carried  out  satisfactorily  using  the  pneumatic  cell  constructed  by  the  author. 
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Synapsis 

Response  function  of  an  optical  system  can  be  obtained  experimentally 
from  the  contrast  of  image  of  a  sinusoidal  grating.  However,  in  practice, 
either  “  finite  grating  ”  with  two  or  three  lines  or  “  variable  grating  ”  whose 
constant  varies  stepwise  from  point  to  point  is  usually  employed. 

The  author  examines  if  the  contrast  of  image  of  such  gratings  gives  true 
response  function  of  the  optical  system  under  test.  In  part  I,  error  of  the 
response  function  obtained  from  “  variable  grating  ”  is  shown  and  the  criterion 
for  the  grating  is  given.  In  part  II,  error  of  the  response  function  obtained 
from  “  finite  greating  ’’  is  shown  analytically. 


1.  Introduction 

In  is  well  known  that  the  response  function  (contrast  transmission  factor) 
of  an  optical  system  can  tie  obtained  experimentally  from  the  contrast  of  the 
image  of  a  sinusoidal  grating  as  follows: 

Transmittance  or  reflectance  of  a  sinusoidal  grating  with  black  background 
across  its  lines  is  given  by 

i4(x)=^l+cosy(jr)) 

where  J{x)  is  a  modulation  function,  its  derivative  being  the  instantaneous  angular 
frequency  uj(x),  that  is: 

%  =“<^>  (>) 

In  a  regular  sinusoidal  grating,  as  is  constant, 

Ax)=aKX 

where  <«o=2x/>io,  >lo  being  the  wavelength. 

If  the  impulse  response  of  an  optical  system,  the  intensity  distribution  of  the 
image  of  an  infinitely  thin  line  at  the  center  of  the  field,  is  h(x),  the  intensity 
distribution  of  the  image  of  the  regular  sinusoidal  grating  is  by  assuming  magnifi¬ 
cation  of  the  optical  system  as  1 


If  we  put 
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h(x)  cos  (M{y-x)dx^ 


5+00  p+oo  p+OO 

h(x)dx=l  and  h{x)cosoM>xdx,  •5®=\  h{x)  sin  (ooxdx  (3) 


B{y) = y[i + Co  cos  otoy + So  sin  owvj 


Accordingly,  contrast  of  the  image  is, 

/?(a>o)=4-"^--"=  ''  Co*(««H3,^)  (  4  ) 

^nux  “T  i3min 

On  the  other  hand,  the  response  function  of  an  optical  system  is  defined  as 
the  Fourier  transform  of  h[x),  that  is. 


HM=  V  Co*{m)+SoKm)  exp  [i  tan'* 


V  VO  C#(aio)J  ^  ^ 

So  the  absolute  value  and  the  phase  of  the  response  function  can  be  obtained  from 
the  contrast  and  the  shift  of  the  image'*. 

But,  in  this  method,  to  find  the  response  function  for  various  m,  numerous 
gratings  of  different  constants  must  be  used.  Moreover,  these  gratings  should  be 
extended  to  infinite  at  both  sides  to  avoid  the  end  effect  (explained  below).  To 
save  these  troubles,  a  grating,  the  constant  of  which  varies  continuously  from 
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point  to  point*  (e.  g.  Fig.  1  (A)),  and  a  grating  of  finite  length*’  (e.  g.  Fig.  1  (B)) 
are  being  used.  But,  with  such  gratings,  we  can  no  longer  put  «>o= const,  in  (4) 
in  the  former  case  and  the  upper  and  lower  limits  of  integral  in  (3)  are  not  infinite 
in  the  latter  case.  Accordingly,  Co  and  So  in  (3)  are  no  longer  equal  to  Co  and 
So  in  (5)  and,  strictly  speaking,  Hiwv)  can  not  be  obtained  from  the  contrast  of 
the  image.  But,  as  such  gratings  are  being  used  extensively  in  practice,  it  is 
very  necessary  to  examine  with  what  degree  of  accuracy  7/(a>o)  is  obtainable  from 
the  contrast  of  the  image  of  these  gratings. 

As  a  variable  grating  corresponds  to  a  frequency  modulated  carrier  and  a  finite 
grating  to  an  amplitude  modulated  carrier  in  terms  of  radio  communication,  we 
shall  examine  the  end  effects  mentioned  above  for  the  case  of  variable  grating  in 
part  I  and  for  the  case  of  finite  grating  in  part  II  by  these  correlations. 

Part  I.  Response  function  obtained  from  the  image  of  a  variable  grating 
§2.  General  formula 

In  the  case  of  variable  grating,  (2)  becomes 

^A(x)cos./(j/-x)dlr]  (6) 

But,  as  h{x)  is  very  small  when  x  is  larger  than  a  certain  value  p,  where  p 
is  about  half  the  width  of  the  image  of  a  thin  line,  we  can  replace  the  upper  and 
lower  limits  of  the  integral  with  p.  Hence,  if  y>p,  then  y>x,  and  j{y—x)  can 
be  expanded  into  Taylor  series  as 

*\y-x)=f{y)-x-f\y)-\-^x*-f'\y) 

If  we  neglect  the  terms  higher  than  x*,  we  get 

J-boe 

A(x)  cos  f{y—x)dx=C  cosf{y)+S  s\x\.f(y) 

-00 

where, 

C=(  h{x)  cos  {f'{y)-x)dx,  S=(  h{,x)  sin  {f'(y)-x)dx 

J  -00  J  -oo 

and  (6)  becomes 

v'C*+^3*cos(/(|/)-v>)]  (7) 

2)  K.  Murata  and  H.  Matsui:  J.  App.  Phys.  Japan  25  (1956)  456. 

3)  O.  H.  Shade:  N.  B.  S.  circular.  526  (1964)  231. 
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where  » 

y=tan-HS/C) 

As  f\x)  is  the  instantaneous  frequency  <w(x)  from  (1),  we  can  see  in  the 
case  of  variable  grating  that  the  contrast  and  the  shift  of  the  image  also  give  the 
absolute  value  and  phase  of  response  function  approximately.  Error  due  to  the 
use  of  variable  grating  is  the  result  of  neglecting  the  terms  higher  than  x*. 

§3.  Modulation  function 

Below,  the  accuracy  of  response  function  obtained  from  (7)  is  examined  with 
some  examples. 


UX1 


Fig.  2.  Modulation  function. 


a  ■*"  IcT  Tf 
C<3 m/yyon 
*■  /To,  o  4  Tfr'm  f 
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In  the  case  of  frequency  varying  as  a  linear  function  of  x  (i.  e.  N.  B.  S,-type 
chart^’ ), 

ai(x)=2mx+l 

The  modulation  function  is  simply  given  from  (1)  as 

flx)=mx^+lx+c  (8) 

But  in  most  variable  gratings  used  in  practice,  frequency  or  pitch  varies  dis- 
continuously  as  is  shown  in  Fig.  1  (A).  Accordingly,  the  modulation  function  is 
not  a  continuous  function  but  a  step  function  as  is  shown  with  full  line  in  Fig.  2. 
At  each  of  the  steps  Pi,Pt,Ps, - Pn,  x  and  Ax)  are  given  as 

n 

x='^ii,  Ax)=2ni:  n:  integer  (9) 

4-1 

As  it  is  hardly  possible  to  treat  such  a  step  function  analytically,  we  shall  use  a 
continuous  function  for  it  as  is  shown  by  a  broken  line  in  the  same  figure  obtained 
by  eliminating  integer  n  as  follows: 

(i)  When  the  pitch  varies  in  geometrical  progression,  we  heve  for  (9) 

j:=xo(l-A-)/(l-/b),  Mx)=2nr. 

where  xo  is  the  first  term  and  k  is  the  common  ratio  of  the  progression. 

The  modulation  function  in  this  case  is,  by  eliminating  n, 

f,{x)=  {2rf\og  k){\og{l-(l-k)xlia)}  (10) 

(ii)  When  the  pitch  varies  in  arithmetical  progression,  we  have  by  putting 

e  the  common  difference  as  d  and  the  first  term  as  xo 

j:=«{2xo+(n— 1)</}/2„  /o(x)=2«r. 

Thus  the  modulation  becomes: 

4(x)={ff(2io-</)/</}{  V  l+&rrf/(2>lo-</)*  -1}  (11) 

In  the  case  ^=1  or  d=0,  both  are  reduced  to 

lim  f,(x) = lim  Mx) = 2ffX/xo 


§4.  Numerical  examples 

The  response  function  is  calculated  in  the  case  of  defocus  by  assuming  h{x) 
to  be  a  square  shape  pulse  of  width  2p,  that  is 

h{x)=\l2p  \x\^p 

hix)=0  \x\>p 

(a)  Geometrical  progression  grating 

The  modulation  function  of  the  grating  in  this  case  is  (10),  but  without  losing 
its  generality,  it  can  be  put  as 


4)  F.  E.  Washer  and  F.  W.  Rosbeny:  J. 

a  aAi,  a 


'liTT  /  /I  -  / 


5*0  — ^ 

ii 


I -A, 


Opt.  Soc.  Am.  41  (1951)  597. 

/I  ■»  nOi  o 

5„  »  ^(jryf  n  Ter<ms 
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fg{x) = A  log  X,  A= 2nl\og  k 

Accordingly,  the  instantaneous  frequency  is 

f,'{x)=Alx=Aly. 

Intensity  distribution  of  the  image  is  by  putting  above  formulas  into  (2) 

B{y)=—^l+^{xlp{l  +  A*)]  cos  (A  log  a;)+/>sin  (A  log  J^)^***} 

So,  the  contrast  is  given  as  follows: 

R(w)=  {sin  /5//>a>}{A/VT+Ai}  V  l+(/to/A)*cot*^ 

where 

,3 = A  log  [(1  +pa)IA)ia -pc»IA)]l2  (13) 

Results  of  calculation  when  ^=1.1  and  ^=0.03,  0.06,'  0.09,  [and  0.12  are 
shown  with  O,  x,  □  and  A,  respectively  in  Fig.  3. 


Rdo) 


Fig.  3.  Contrast  of  the  image  of  sintiisoidal  variable  grating  (geometrical  progression). 

Full  line  denotes  correct  response  function. 

O,  X,  □  and  A  show  calculated  contrasts  for />= 0.03, />= 0.06, />= 0.09  and  ^=0.12, 
respectively. 

The  response  function  in  this  case  should  be 

f/(a>o)=sin  pmlpon  (14) 

The  Full  line  in  Fig.  3  shows  this  correct  response  function. 

If  {pa)IAY  in  (13)  is  assumed  to  be  negligibly  small,  we  can  put  approx¬ 

imately,  and  we  get 


/?(<«)= sin  pwlpci). 
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Thus  the  contrast  of  the  image  of  a  variable  grating  becomes  the  correct  response 
function. 

Calculated  values  of  the  error 

(15) 

for  various  k  are  given  in  Fig.  4  which  affirms  that,  for  measuring  Him)  with 
the  accuracy  of  5%  when  Pm<!2.Tt,  a  grating  of  *<1. 15  should  be  used.  Error 
in  phase  of  the  response  function  is  negligibly  small  in  this  case. 


Fig.  4.  Percentage  error  of  response  function  obtained  from  contrast  of  the  image*of 
sinusoidal  variable  grating  (geometrical  progression). 

(b)  Arithmetical  progression  grating 

The  Modulation  function  of  the  grating  is  given  by  (11)  which  can  however 
be  put  in  general  as 


f.(x)=y/xip. 
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Intensity  distribution  of  the  image  in  this  case  is  given  from  (2)  as 
B{y)={\+{^lp)[g{y+p)-g(,y-p)\)l2 

where 


{/(/)=cos^+^sin6>, 

Results  of  calculation  of  the  contrast  therefrom,  when  Vl/’^=10;r  and  />=0. 1,  0.2, 
and  0.3,  are  shown  in  Fig.  5  with  O,  X  and  □,  resiiectively.  In  this  figure  the 
full  line  denotes  the  correct  value  of  the  response  function  (14).  Comparing  with 


R(a)) 


Fig.  5.  Contrast  of  the  image  of  sinusoidal  variable  grating  (arithmetical  progression). 
Full  line  denotes  correct  response  function. 

O.  x»  and  □  show  calculated  contrasts  for  p=0.1,  p=0.2  and  />=0.3,  respectively. 


Fig.  3,  we  can  see  that  the  error  in  this  case  is  a  little  larger  than  that  with  the 
geometrical  progression  grating  when  xo  (initial  wave  length)  is  the  same  in  both 
gratings.  Calculated  values  of  the  percentage  error  by  (15)  for  various  d  are  given 
in  Fig.  6  in  which  the  error  of  amplitude  is  shown  with  full  lines  and  that  of 
phase  with  broken  lines.  It  is  seen  that  the  error  depends  upon  the  width  p,  and 
to  get  the  response  function  as  accurate  as  in  the  case  of  geometrical  progression 
grating,  a  grating  of  <///><0. 1  should  be  used. 

(c)  N.  B.  S.-type  grating 

In  this  case  the  modulation  function  is  given  by  (8),  the  integral  in  (6)  is 
therefore  a  Fresnel  integral  and  it  is  hardly  possible  to  find  C(<u)  and  S(a>)  analy¬ 
tically.  But,  if  we  assume  p<y,  it  becomes 


L 
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Fig.  6.  Percentage  error  of  response  function  obtained  from  contrast  of  the 
image  of  sinusoidal  variable  grating  (arithmetical  progression).  Full  lines 
show  error  for  absolute  value.  Broken  lines  show  error  for  phase. 


R(w)=  {sin  2//>^o>T+cot  A«'/<'4** 

Accordingly,  percentage  error  of  the  response  function  is  approximately 

^=P*m^l2, 

except  when  pw=nn,  where  n  is  an  integer.  As  this  expression  of  the  error  does 
not  contain  w,  the  response  function  with  constant  percentage  error  up  to  a  very 
high  frequency  can  be  evaluated  which  is  a  grate  merit  of  N.  B.  S.  -type  grating. 
In  Fig.  7,  results  of  calculation  of  R(a»)  are  shown  when  m=ff/l(X).  The  correct 
value  of  the  response  function  is  shown  by  the  full  line  on  which  x ,  O  denote 
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X  :  p  .  2.8 
a:  p.3.S 
o  :  p>4.2 


Fig.  7.  Contrast  of  the  image  of  sinusoidal  N.  B.  S.-type  variable  grating.  Full  line 
denotes  correct  response  function. 

X,  A  and  O  show  calculated  contrasts  for  p=2.8,  p=3.5  and  p=4.2,  respectively. 

the  calculated  values  p=2.8,  3.5  and  4.2  respectively.  As  is  expected,  the  differ¬ 
ences  between  the  values,  calculated  and  correct,  are  kept  very  small  up  to  a 
very  high  frequency. 

(d)  Square  wave  variable  grating 

The  result  of  the  examination  for  the  case  of  variable  square  wave  grating 
shows,  in  general,  that  the  error  is  only  a  little  larger  than  that  with  the 
sinusoidal  grating  shown  above.  If  the  response  function  obtained  from  the  square 
wave  and  the  sinusoidal  grating  are  denoted  as  rim)  and  Rim)  respectively,  the 
Coltman  formula®’  gives  their  relation  as 

•  •  •] 

It  became  clear  now  that  the  errors  in  measurement  of  Rim)  and  rim)  are  of  the 
same  order  for  the  same  m.  Hence  the  calculated  value  of  Rim)  from  the  experi¬ 
mental  results  of  rim)  by  the  above  formula  is  less  accurate  than  that  obtained 
direct  from  the  sinusoidal  grating. 

Therefore,  a  square  wave  grating  with  a  constant  which  is  smaller  than  that 

t 

of  a  sinusoidal  grating  has  to  be  used. 

Part  II.  Response  function  obtained  from  the  image  of  finite  grating 
§5.  General  formula 

Intensity  distribution  of  a  finite  grating  as  is  shown  in  figure  1(B)  is  given 
in  general  as 

A  ix) = kix)  -h  g  ix)  cos  oHiX 

where  kix)  is  a  constant  term  and  gix)  the  amplitude  modulation  function  given  as 
5)  J.  W.  Coltman:  J.  Opt.  Soc  Am.  44  (1954)  468. 
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Ar(x)=a+ft=(H-6)/2,  g{x)=a  \x\^W 

k{x)=c,  g(x)=0  \x\>W 

where  2W  is  the  total  length  of  the  grating,  a  the  amplitude,  b  and  c  the  in¬ 
tensities  of  the  background  and  both  sides  of  the  grating  which  is  called  the 
surround  by  W.  H.  Steel*’,  respectively  (Fig.  8).  a,  b  and  c  are  normalised  as 
2a-\-b=\  and  c:^\. 


MX) 


The  intensity  distribution  of  the  image  is  given  from  (2)  as 

B{y)=V^ k{y-x)h(x)dx+[*  g{y-x)h(x) cos <jM,(y-x)dx=Bi{y)+B2(y)  (16) 

J-oo  J-oo 

The  first  term  is  the  image  of  a  single  slit  with  half  width  W,  and  can  be 
calculated  when  the  impulse  response  is  given. 

The  intensities  at  the  points  y=0  and  y=yo  in  the  image  plane  are  given  as 
shown  in  Fig.  9(B)  as 


(i)  y=0,  W^p,  Bi=c\  h{x)dx+(\  h{x)dx=c+(=(l+b)l2 

J-p  J-p 

w<p,  Bi=c+(a-do) 


where 


^0  = 


■-IT  pp 

h{x)dx+\  h(x)dx 

-p  JW 


(17) 


6>  W.  H.  Steel;  Optica  Acta.  3  (1956)  456. 
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(B)  Image  of  single  slit 
Fig.  9.  Intensity  distribution  of  single  slit. 


0 


(ii)  If =2/0  H^^/>+l/o 

W^O+i/o 


fit=c+€=(l+6)/2 
.6i=c+f(l— tf) 


(•-•O'  fKO+P 

where  6=\  h{x)dx+\  h{,x)dx 
J»o-p  J'*' 

As  to  the  second  term  of  (16), 'if  we  put  the  Fourier  transforms  of  h{x)  and 
g{,x)  as  //(«»)  and  G{u),  then 


where  Re  represents  the  real  part  of  complex  quantities. 

If  W  is  assumed  sufficiently  large,  u  will  be  limited  in  a  small  domain  2a, 
and  //(«>o+m)  can  be  expanded  into  Taylor  series  as 


and  the  amplitude  of  Bt  becomes  as 
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G{u)e^'‘*du  +  G(m) •  u  •  •  •  •  •  (18) 

where  H'iwt),  //"(a*)*  • -are  the  successive  derivatives  of  H{w9).  But,  if  H\uh>) 
and  a  are  not  so  large,  some  preliminaly  calculations  showed  that  the  second  term 
has  little  effect  on  the  result  and  can  be  neglected  leaving  only  the  first  term  to 
be  considered. 

The  first  term  of  (18)  is  well  known  in  terms  of  radio  communication  as  the 
wave  form  of  output  signal  of  square  wave  input  transmitted  through  an  ideal 
band-pass  filter  with  width  2a,  and  we  get’’ 

-^//(ai.)[Sf(a(y  -f  ly  ))-S»(«(y-  ))]= a//(a>oXl  -  jj)  (19) 

where  Si(i)=( -  dz, 

Jo  2 

and  the  intensities  at  i/=0  and  l/=i/o  are  from  (17)  and  (19) 
l/=0:  — ^)+fl//(a>oKl — >?•) 

y=yo:  B=c+f(l-3)+o//(€«*o)(l-iy)e‘*»»« 

where  r«=l— -^1(01^). 

Hence,  the  contrast  of  the  image  obtained  from  the  intensities  at  i/=0  and 
yt=7:lw  is  approximately 

As  is  seen  from  this  formula,  although  Riwtt)  does  not  give  the  correct  response 
function  (due  to  what  we  shall  call  the  “end  effect”),  the  error  JH{wo)=R{m)— 
Him)  is  the  smaller  the  larger  so  a  grating  of  a  large  number  of  lines  should 
be  used  to  get  a  better  result. 

The  error  also  becomes  small  by  selecting  a  suitable  value  of  <.  This  is, the 
reason  why  a  grating  with  gray  surround  is  more  effective  to  reduce  the  “end 
effect”  than  that  with  black  surround®’. 

If  the  attenuation  of  response  function  Him)  is  montonous  as  in  the  case  of 
an  aberrationless  ideal  optical  system,  the  grating  with  gray  surround  of  c=0 
(that  is,  if  6=0,  reflection  of  transmission  of  the  surround  is  50%)  gives  the  best 
result,  because  in  this  case  and  1^  are  very  anall  as  compared  with  dt  and  8. 
However,  in  the  case  of  a  practical  optical  system  with  aberrations,  >;•  and  1;  are 
not  negligibly  small;  the  “end  effect”  can  be  reduced  by  selecting  a  suitable 
value  of  (  as  shown  by  the  following  example. 

7)  e.  g.  S.  Goldman:  Frequency  Analysis,  Modulation  and  Noise.  (McGraw-Hill,  N.  Y. 
(1948)  72. 


§6.  Numerical  example 


As  an  example  of  the  “end  effect”  given  by  (20),  the  case  in  which  the 
impulse  respone  is  of  a  square  shape  as  is  given  by  (12)  is  shown.  The  contrasts 
of  the  finite  sinusoidal  grating  are  calculated  for  the  case  in  which  the  number 
of  lines  are  N=2,  3,  and  4,  and  c=0,  0.3  and  0.5,  with  a =0.5,  6=0  so  f=0. 5— c 
and  /to: =0.35.  The  results  are  plotted  by  chain  lines  for  c=0  (grating  with  black 
surround),  by  dotted  lines  for  c=0. 3  and  by  broken  lines  for  c=0. 5  in  Fig.  10. 


Fig.  10.  Contrast  of  the  image  of  sinusoidal  finite  grating  with  2,  3  and  4  lines  for 
various  contrasts  of  surround.  Full  line  denotes  correct  response  function,  chain 
lines  denote  calculated  contrast  for  c=0,  dotted  lines  for  c=0.3,  broken  lines  for 
c=0.5 

The  correct  response  function  qf  this  case  is  given  from  (14)  and  plotted  by 
the  full  line  in  the  same  figure. 

If  the  grating  of  c=0. 3  is  used,  R{m)  becomes  nearer  to  the  correct  value 
ff(a)o)  than  using  either  black  surround  (c=0)  or  gray  surround  of  c=0. 5.  But, 
it  must  be  remembered  that  the  value  of  e  which  gives  the  best  result  depends 
upon  the  type  of  the  impulse  response. 

However,  we  can  see  from  this  figure  that  the  results  obtained  from  the  con¬ 
trast  of  the  image  are  correct  in  all  the  cases  up  to  about  (N— l)th  zero  of  the 
response  function. 

Intensity  distributions  of  the  image  of  black  and  gray  surround  gratings 
(c=0. 5)  when  /V=3  are  shown  in  Fig.  11  (A)  and  (B)  for  various  values  of  pwo. 
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Fig.  11.  Intensity  distribution  of  the  image  of  sinusoidal  finite  grating  with  3 
lines  for  various  width  of  impulse  response. 

(A)  in  the  case  of  black  surround  (c=0). 

(B)  in  the  case  of  gray  surround  (c=0.5). 

The  results  of  calculation  of  the  contrast  using  the  values  at  a}oy=0  and  oJoj/=r 
in  this  figure  are  plotted  with  O  and  x  in  Fig.  10.  They  are  in  good  agreement 
with  the  calculated  value  from  (20). 

Conclusion 

Examination  was  made  if  the  contrast  of  the  image  of  “variable  grating” 
and  “finite  grating”  gives  the  correct  response  function.  As  to  the  case  of 
variable  grating,  it  is  shown  generally  that,  when  the  instantaneous  frequency 

=a>(ar)  is  taken  as  spatial  frequency,  the  contrast  of  image  rejwesents  approxi- 
dx 

mately  the  correct  response  function.  This  is  confirmed  in  Part  I  by  numerical 
calculations  with  three  kinds  of  sinusoidal  gratings,  that  is,  the  cases  when  the 
pitch  varies  in  geometrical  and  arithmetical  progressions  and  the  case  of  N.  B.  S. 
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type  grating.  Their  modulation  functions  and  constants  of  grating  with  which 
the  error  is  less  than  \%  up  to  4th  zero  of  the  response  function  for  a  square 
shape  impulse,  etc.  are  tabulated  in  Table  1. 


Table  1.  Sinusoidal  variable  grating 


Type  of  varying 
constant 

Modulation 

function 

Upper  limit 
of  constant 

N.  B. 

Geometrical 

1  progression 

A  logx 

*<1.05 

Error  is  inde¬ 
pendent  width  of 
impulse 

j  Arithmetical 
progression 

dip  <0.02 

j 

N.  B.  S. 
type 

ntj^+lx 

mpi‘-<0.\A 

Error  is  inde¬ 
pendent  of  spatial 
frequency 

NOTE:  The  upper  limit  of  constant  represents  the  constant  with  which 
error  of  response  function  is  less  than  \%  up  to  4th  zero  of  response 
function  for  a  square  shape  impulse.  This  limit  of  spatial  frequency 
corresponds  to  the  limit  of  5  line  finite  grating. 

The  error  with  a  square  wave  grating  is  of  the  same  Oder  as  that  with  a 
sinusoidal  grating  of  the  same  constant.  Therefore,  in  oder  that  the  accuracy,  of 
the  response  function  obtained  from  a  square  wave  grating  becomes  equal  to  that 
of  a  sinusoidal  grating,  the  constant  of  the  square  wave  grating  should  be  much 
smaller  than  that  of  the  sinusoidal  grating.  / 

In  part  II,  the  error  with  a  finite  grating  is  considered  analytically.  The  re-  * 

suits  show  that  the  error  depends  upon  the  number  of  lines  and  the  contrast  of 
the  grating.  If  the  grating  has  N  lines,  the  correct  response  function  will  be 
obtained  up  to  the  (iV— l)th  zero  of  the  response  function  for  a  square  shape  im¬ 
pulse.  It  is  also  shown  that  the  gray  surround  is  more  profitable  than  the 
black  surround. 

When  the  so  called  “  Siemens’  Star  ”  is  used  as  the  chart  for  measuring  the 
response  function,  there  also  occur  errors  of  the  same  kind  due  to  the  continuous 
change  of  the  grating  constant  along  the  radius  of  the  star.  As  to  the  error  in 
this  case  we  have  an  excellent  computation  by  K.  P.  Miyake*’. 
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Abatrsct 

Arc  spectra  in  air,  nitrogen,  oxygen  and  argon  using  pure  aluminum  of 
high  grade  as  two  electrodes  are  photographed  by  a  quartz  and  a  glass  prism 
spectrographs.  Many  new  bands  and  atomic  lines  appear  between  2500  A  and 
8800  A.  Comparison  of  spectra  in  various  gases  suggests  emitters  of  new  bands. 

All  the  heads  of  {A*I—X*S)  system  of  aluminum  monoxide  reported  by  Roy 
are  examined  and  about  20  new  heads  are  added  to  this  system  extending  the 
vibrational  levels  up  to  i;''=12  of  X*I  and  v^=16  of  A^I  state.  It  is  indicated 
that  r"=9  level  of  X*S  is  raised  about  10  cm-*  on  account  of  perturbation  by 
a  certain  state. 

In  the  near  infrared  spectrum  several  atomic  lines  of  A1  i  are  analyzed 
and  the  height  of  n^F  (n=8,  9  and  10)  levels  are  estimated. 

1.  Introduction 

Among  aluminum  compounds  excepting  hydride  and  halide  only  aluminum 
monoxide  is  ever  known  to  display  its  spectrum,  and  only  two  electronic  states, 
A*I*  and  have  so  far  been  confirmed.  Since  1866  there  have  been  many 

reports*' ■*’  on  the  (A*!"— X*!")  system,  which  appears  in  blue-green  region  fon^g 
clear  groups  of  heads,  and  in  1925*’  their  emitter  was  proved  to  be  the  diatomic 
oxide  molecule,  AlO.  In  1939  D.  Roy*’  made  precise  measurements  of  the  heads 
of  this  system  including  many  new  ones  between  4150-57(X)A.  As  the  result  of 
his  analysis  the  vibrational  levels  of  and  state  have  been  extended  re¬ 
spectively  to  t>'=14  and  »"=11.  . 

Later  the  predissociation  in  the  (0—0)  band  of  this  system  at  A=129  was 
observed  by  C^heur-Dehalu^,  which  was  supported  by  B.  Rosen*’.  Rosen  proposed 
the  predissociation  in  the  lower  *£  state  on  account  of  perturbation  by  a  certain 
repulsive  state,  and  concluded  that  all  the  bands  with  v"^9  are  predissociated. 

1)  R.  Mecke:  Phys.  Zeit.  26  (1925)  217. 

2)  G.  Eriksson  and  E.  Hult^:  Z.  Phys.  34  (1925)  775. 

3)  W.  C.  Pomeroy:  Phys.  Rev.  29  (1927)  59. 

4)  M.  K.  Sen:  Indian  J.  Phys.  11  (1937)  251. 

5)  D.  Roy:  Indian  J.  Phys.  13  (1939)  231. 

6)  R.  Mulliken:  Phys.  Rev.  26  (1925)  561. 

7)  F.  Coheur-Dehalu:  Bull.  Acad.  Roy.  Belgium  23  (1937)  604. 

8)  B.  Rosen:  Phys.  Rev.  68  (1945)  124. 
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From  this  proposition  it  follows  that  one  of  the  two  atoms  resulting  from  the 
dissociation  of  AlO  in  X*I  state  must  be  in  excited  state;  the  products  of  the 
dissociation  are  probably  A1(*F*)  and  0{^D).  And  according  to  this  view,  the  dis¬ 
sociation  energy  should  be  as  low  as  0. 93eV.  Recently,  however,  Lagerqvist, 
Nilssson  and  Barrow*',  after  examining  the  (0-0)  band  of  the  system 

with  a  high  resolving  power  spectrograph,  stated  that  the  observed  predissociation 
at  P(129)  is  really  an  instrumental  effect  and  that  all  the  other  reported  predis¬ 
sociations  in  this  system  might  have  been  illusory. 

As  regards  other  systems  of  aluminum  monoxide,  F.  P.  Coheur  and  B.  Rosen 
published  a  report  in  Using  the  exploding  wire  technique,  they  found 

many  heads  in  the  region  of  2500—  3500  A  and  tried  to  analyze  them.  In  their 
result  there  is  (B—X^I)  system**'  and  a  reasonable  Deslandres  table  with  28  heads 
was  formed.  However  it  seems  desirable  to  examine  this  system  further  because 
the  heights  of  vibrational  levels  of  lower  state  with  v"^5  are  lower  than  those 
obtained  from  (A^I—X^I)  system. 

It  is  also  interesting  to  search  for  a  *//  state  corresponding  to  A^fh  state  of 
BO  and  for  this  purpose  an  investigation  in  the  near  infra-red  seems  necessary. 
Moreover  the  spectra  of  other  aluminum  compounds  such  as  AIN  and  AIS  qre 
interesting  from  the  point  of  view  of  molecular  structure.  It  seems  also  worthwhile 
searching  for  the  spectrum  of  AU  because  that  of  Bz  has  already  been  investigated**’ . 

Spurred  by  this  interest  the  author  sought  the  methods  of  exciting  aluminum 
compounds.  The  arc  between  pure  aluminum  electrodes  in  various  gases  appears 
fairly  useful  for  this  purpose.  The  arc  in  air  in  this  method  displayed  the  spectra 
of  many  bands  in  the  region  of  25(X)— 35()0A,  and  {A^I—X^I)  system  of  AlO  was 
extended  to  4050  A  revealing  one  more  sequence  than  has  ever  been  reported. 
The  arc  spectra  in  nitrogen,  oxygen  and  argon  with  pure  aluminum  electrodes 
were  photographed,  and  some  of  the  new  heads  and  atomic  lines  appearing  in 
them  were  analyzed. 

2.  Experiments 

The  arc  was  run  from  D.  C.  100  V  supply  through  a  resistance  of  20  ohms. 
Pure  aluminum  of  high  grade  was  employed  as  the  two  electrodes  in  order  to 
avoid  the  spectrum  of  any  impurity.  The  aluminum  was  graded  99.994%  pure 
with  Fe,  Cu  and  Si  as  impurities,  each  of  which  being  less  than  0. 002%  with  just 

9)  A.  Lagerqvist,  N.  E.  Lennart  Nilsson  and  R.  F.  Barrow:  Proc.  Phys.  Soc.  69  (1956)  356. 

10)  F.  P.  Coheur  and  B.  Rosen:  Mem.  Soc.  Roy.  Sci.  Li^e  (1941)  405. 

11)  G.  Herzberg:  Spectra  of  Diatomic  Moleules  (1950)  p.  504  (Van  Nostrant  Comp.,  New 
York). 

12)  A.  E.  Douglas  and  G.  Herzberg:  Canad.  J.  Res.  A  18  (1940)  165. 
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a  trace  of  Mn.  Several  shapes  of  electrodes  were  tried  in  search  for  the  best  to 
prolong  the  arcing  and  at  the  same  time  to  facilitate  the  display  of  spectra  of 
many  bands.  A  rod  of  10  mm  in  diameter  as  the  lower  negative  electrode  and  a 
rolled  foil  of  0. 1  mm  in  thickness  or  two  or  three  wires  of  1  mm  in  diameter 
twisted  together  as  the  upper  electrode  were  found  fairly  satisfactory  for  the 
purpose.  It  was  necessary  to  keep  the  electrodes  free  from  oil  or  fat,  for  aluminum 
seemed  to  adsorb  them  very  easily  on  its  surface  causing  strong  CN  bands  to 
appear  impeding  the  observation.  The  discharge  tube  was  made  from  two  liter 
flask  and  one  of  the  electrodes  was  operated  by  a  bellows. 

The  pressure  of  gas  was  kept  at  almost  760mmHg,  for  in  low  pressures,  gas 
appeared  easily  excited  to  show  its  spectrum.  The  employed  gases  were  nitrogen 
and  argon  of  high  grade  and  commercial  oxygen  in  bomb.  Nitrogen  was  further 
purified  by  passing  it  through  heated  copper  flinders  to  remove  even  a  trace  of 
oxygen  to  avoid  the  other  spectra  being  obscured  by  the  blue-green  system  of  AlO. 
Arc  in  purified  nitrogen  or  argon  never  goes  out  provided  the  distance  between 
the  electrodes  is  properly  kept.  But,  since  arc  in  oxygen  dies  out  in  no  time, 
photographing  of  the  spectrum  of  AlO  was  made  with  arc  in  air.  The  following 
seemed  worth  remembering  for  the  analysis  of  the  result:  both  electrodes  and 
walls  of  the  discharge  tube  were  whitened  with  aluminum  oxide  (AhOs)  powder 
when  even  a  trace  of  oxygen  existed  in  the  discharge  tube;  they  were  coloured 
grayish  black  with  aluminum  nitride  (AIN)  powder  when  purified  nitrogen  was 
used;  they  remained  clean  in  argon. 

The  spectra  in  visible  and  near  infra-red  regions  were  photographed  by  a 
glass  spectrograph  composed  of  three  60°  prisms  of  6  cm  base  and  a  camera  lens 
of  50  cm  focal  length.  The  dispersion  was  about  46  A/mm  at  6(X)0A.  For  ultra¬ 
violet  region  a  Hilger  Ei  quartz  spectrograph  was  used.  The  photographic  plates 
employed  were  Fuji  process  and  Neopan  Type  for  ultra-violet  and  visible  regions 
respectively,  Kodak  I-N  for  near  infra-red  region. 

Iron  arc  lines  were  taken  as  reference  standards  through  the  whole  region. 
The  probable  error  in  wavelength  determination  was  about  ±0. 1 A  in  the  range 
observable  with  Ei  spectrograph,  ±0.3  A  in  the  visible  region.  But  it  amounted 
to  ±1.2  A  in  the  near  infra-red  region. 

3.  Results  and  Discussion 

(i)  Visible  and  near  ultra-violet  region 

This  is  the  region  where  {A*1—X^S)  system  of  AlO  dominates.  The  arc 
spectrum  in  air  showed  almost  twenty  new  heads.  Removal  of  disturbing  Fe  atomic 
lines  and  heads  of  CN  violet  system  enabled  a  new  sequence  of  Jt;=5  between 


Table  1.  Deslandres  table  of  system  of  AlO 
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The  values  in  Gothic  are  new,  and  those  marked  with  *  are  revised. 
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865.8 

4^.20 

(5) 

21614.6 

not  clear 
in  our 
spectrum 

5016.9 

(0) 

19927.1 

5464.0  1 

(4d)! 

18296.5; 

1  79^41  i  788.2 

788.9 

i  4139.6  ,  1  4296.44 

(1)1  !  (3) 

1  24156.ll  881.51  23268.6 

865.1 

4462.33 

(4) 

22403.5 

861.0* 

4640.7* 

(4d) 

2154Z5 

. 

5024.4 

(1) 

19897.3 

* 

1  1 

783.8* 

782.5 

in  the 
Aadc  of 

(4-0) 

4311.5* 
1  (3) 

23187.3 

863.1* 

4478.2* 

(4d) 

22324.2 

1(4650.30) 

(826.2)  (21498.0 

i 

5032.6 

(0) 

19864.91 

i 

:  774.7 

(765.6 

1 

[ 

4172.1 

!  (2) 

1  23962.C 

i 

oTerlnp 

with 

(4-1) 

,  4490.38 
(3) 

j  22263.e 

i 

823.2 

1  4662.8 

(1) 

t|  21440.2 

1 

766.1 

1  764.2 

1 

i  4196.8 
(2) 
2S855.1 

825.4 

4341.0 

(2) 

28029.' 

825.1 

1  4502.3^ 

i  (2 

S  22204.; 

■ 

y 

1 

t 

1 

r  ' 

i 

in  the 
!ahade  o 

1  (0-2) 

in  the 
shade  o 

1  (5-2) 

f; 

4512.3 

(1) 

!  ffil55.5 

1 

! 

j 

750.2! 

1 

! 

1 

1 

1 

1 

4216.5 

(2 

1  23709. 

)‘ 

6.  803. 

;  4364.5  ' 

(2d)! 
9  22905.7 
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4041 A  and  4217  A  to  appear.  In  this  sequence  the  bands  with  high  vibrational 
quantum  numbers  have  fairly  strong  intensities  and  the  band  of  (16—11)  is  clearly 
observed.  The  result  of  analysis  is  shown  in  Table  1,  in  which  the  newly  found 
heads  are  in  Gothic  and  numbers  in  brackets  are  the  estimated  intensities.  All 
the  heads  reported  by  D.  Roy*’  were  examined,  and  the  wavelengths  of  those 
heads  with  which  the  present  result  agree  within  experimental  error  are  given  in 
Roy’s  figures  without  any  sign,  for  his  measurements  appeared  more  precise  than 
the  author’s  on  account  of  better  spectrograph.  The  values  which  could  not  be 
exmined  are  in  parentheses.  At  the  wavelengths,  5445. 05  A,  5216. 20  A  and 
4476. 02  A,  which  were  reported  by  Roy*’  as  the  positions  of  heads  of  (7—9),  (8—9) 
and  (12—9)  bands  respectively,  no  such  indication  was  found.  However  about  2  A 
on  the  longer  side  of  each  of  them  there  are  what  are  to  be  regarded  as  heads.  The 
average  difference  of  wavenumbers  between  Roy’s  and  the  author’s  is  about  10  cm"*, 
which  far  exceeds  the  experimental  error.  Similarly  what  exists  at  4638. 77  A  that 
is  reported  by  him  as  the  head  of  (11—9)  band  looks  hardly  acceptable  as  a  head, 
and  what  exists  slightly  (1. 9  A)  on  the  longer  side  of  it  appears  more  like  a  head. 
Therefore  the  revised  values  are  listed  in  Deslandres  table  with  asterisks.  Taking 
acount  of  all  the  revision  and  the  positions  of  new  bands  with  v"=9,  it  is  con¬ 
cluded  that  the  vibrational  level  of  v"=9  of  X^I  state  is  perturbed  and  raised 
about  10  cm"*  by  a  certain  state.  All  the  heads  of  bands  with  v"=9  seem  a  little 
more  diffuse  than  the  others.  The  analysis  of  the  vibrational  levels  was  extended 
up  to  v"=l2  of  state,  but  no  appreciable  perturbation  to  the  levels  of  t;"=10, 
11  and  12  was  found,  for  the  calculated  wavenumbers  using  Roy’s  foamula  agree 
with  those  observed  within  the  experimental  error.  It  should  be  noted  that  the 
head  of  (12—8)  band  does  not  appear  at  the  reported  position*’,  and  the  revised 
value  is  listed  also  with  an  asterisk. 

At  the  tail  of  Jv=0  sequence,  some  heads  seem  to  overlap  some  rotational 
lines,  and  the  positions  of  abnormaly  strong  lines  coincide  with  the  calculated 
values  of  heads  of  (8-8),  (9-9),  (10-10),  (11-11)  and  (12-12)  bands.  The  ap- 
perance  of  these  bands  seems  reasonable,  for  they  might  probably  represent  sub- 
sidary  intensity  maxima  on  Deslandres  table. 

The  arc  in  nitrogen  or  in  argon  shows  no  band  spectrum  in  this  region. 
Atomic  lines  of  aluminum  are  easily  excited  in  the  arc  in  nitrogen,  and  several 
new  lines  appear  when  pure  aluminum  of  high  grade  is  employed  as  electrodes. 
Although  they  are  not  capable  of  being  combined  with  the  known  atomic  energy 
levels  of  A1 1**’ ,  they  might  possibly  be  the  atomic  lines  of  A1 1,  since  no  line  of 

13)  C.  E.  Moore:  Atomic  Energy  Levels  I,  (1949)  p.  124-143  (Circular  of  National  Bureau 
,of  Standards  467,  Washington). 


108 


Midori  Shimauchi 


A1  or  A1  appeared  over  the  whole  region  investigated.  For 

reference,  some  values  of  them  are  listed  in  Table  2,  and  their  photographs  are 
shown  in  Fig.  1,  in  which  the  arc  spectra  in  air  and  in  nitrogen  are  placed  side 
by  side. 


Table  2.  Unanalyzed  lines  of  A1 


Int. 

V 

•  Reference 

* 

^20.5 

2 

24865.5 

* 

4023,7 

2 

24845.7 

4207.9 

1 

23758.1 

* 

5073.9 

4d 

19703.2 

5108.0 

3d 

19571.7 

slightly  shades 
towards  red 

5173.2 

3 

19325.0 

5184.4 

4 

19283.3 

5588.7 

1 

17888.3 

♦ 

5594.1 

0 

17871.0 

♦ 

6122.6 

1 

16328.4 

6162.6 

2 

16222.4 

6440.8 

1 

15521.7 

6464.4 

1 

15465.1 

All  the  lines  listed  are  observed  most  clearly  in  the  arc  spectrum  in 
nitrogen,  while  those  marked  with  *  appear  sometimes  in  the  arc  in  air. 


(2)  Near  infra-red  region 

In  this  region,  the  spectrum  in  air  or  in  oxygen  is  obscured  by  a  continuous 
spectrum,  whereas  that  in  nitrogen  is  most  interesting.  In  this  spectrum  no  atomic 
line  of  nitrogen  appears,  while  in  the  arc  in  argon  atomic  lines  of  argon  are 
easily  excited.  In  Fig.  2  the  arc  spectra  in  nitrogen  and  argon  are  shown  side  by  side 
for  the  convienience  of  comparison.  In  the  arc  in  nitrogen  the  diffuse  (4*P— n*Z)) 
and  the  sharp  {A^D—n^S)  series,  which  have  been  observed  only  in  some  specific 
discharges^’ appeared  clearly  vrith  more  intensity  than  in  the  discharge. 
Three  new  lines  listed  in  Table  3  are  readily  analyzed  as  the  extension  of  these 
series.  Their  wavelengths  estimated  from  “  Atomic  Energy  Levels'”  ”  are  in 
accord  with  those  observed.  The  lines  of  fundamental  (3*D— «*F)  series  appear 
strongly  with  shade  about  forty  to  seventy  cm“'  towards  red  in  the  arc  spectrum 
in  nitrogen,  while  they  have  no  shade  when  the  atmosphare  of  arc  is  argon. 
This  series  was  extended  up  to  n=10.  The  observed  wavelengths  of  lines  and 

14)  F.  Paschen  and  R.  Ritschl:  Ann.  Phys.  [Leibzig]  (5)  18  (1933)  886. 

15)  G.  R.  Harrison:  Af.  1.  T.  Wavelength  Tables  (1939)  (John  Willey  and  Sons,  Inc.,  New 
York) 

16)  Spectraltabellcn  (1955)  (Veb  Verlag  Technik,  Berlin). 
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(a)  A1  arc  in  argon 

(b)  A1  arc  in  nitrogen 


Fig.  2 


Table  3.  Analyzed  lines  of  A1  i 


1  ^ob* 

1 

Int 

^cal 

Veal 

Transition  j 

!  7605.6 

3 

7605.3 

13144.06 

(4*Ps/,-8*Di/,) 

!  7563.4 

2 

7562.9 

13218.83 

(4*P3/*-9»Sv,) 

7555.1 

1 

7553.8 

13234.66 

(42Pi/,-9*Si/,) 

1 

Table  4.  Rydberg  Series  of  A1  i 


n 

iobi 

Int. 

Votl 

height  of 
ri^F 

quantum 
defect,  fi 

4 

41318.74* 

0.029 

5 

(8772.88) 

8773.1 

10 

(11396.63) 

11395.4 

43831.08* 

0.033 

6 

(7835.33) 

7835.3 

9 

(12759.20) 

12759.2 

45194.65* 

0.035 

7 

(7361.59) 

7362.8 

8 

(13580.27) 

13578.9 

46015.73* 

0.037 

8 

7085.2 

5 

14110.4 

46545.85 

0.043 

9 

6907.3 

3 

14473.5 

46908.95 

0.051 

10 

6785.3 

_ 

14733.7 

47169.15 

0.056 

The  values  in  brackets  are  taken  from  the  report  of  Paschen  and  RitschP^’,  and 
those  marked  with  ♦  from  “Atomic  Energy  Levels”**. 
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estimated  heights  of  «*F  levels  are  listed 
in  Table  4.  Quantum  defect,  (i,  was  cal¬ 
culated  using  Rydberg  constant  of  aluminum, 
109735  cm**,  which  was  obtained  from  the 
following  formula: 


A :  atomic  weight 
Af,:  mass  of  the  proton 
w:  mass  of  the  electron 
/?oo:  Rydberg  constant  (infinite  mass). 
The  variation  of  ft  with  the  increase  in  n 
is  plotted  in  Fig.  3.  The  fi  values  of  newly 
estimated  states,  S^F,  9*F  and  10*F,  are 
precisely  on  the  fextension  of  the  curve 
which  is  determined  by  the  points  calculated 
from  the  heights  of  already  known  states,  4*F,  5*F,  6*F  and  7*F‘*’ . 

There  are  several  new  heads  of  bands,  whose  shades,  wavelengths  and  wave- 
numbers  are,  for  reference,  listed  in  Table  5.  It  is  possible  that  these  bands  are 
emitted  from  aluminum  nitride  or  aluminum  molecule,  since  they  ap(>ear  clearly 
when  purified  nitrogen  is  used  as  the  ambient  gas,  while  even  the  strongest  head 
of  {A*S—X*I)  system  of  AlO  does  not  appear.  The  arc  in  argon  reveals  no  bands. 

Some  of  the  heads  observed  between  2500  A  and  3500  A  coincide  with  those 
reported  by  F.  P.  Coheur  and  B.  Rosen**’ .  They  are  certainly  emitted  from  oxide 
of  aluminum,  because  they  appear  in  the  arc  spectrum  in  air  or  in  oxygen  and 


Fig.  3.  Quantum  defect  of  tfiF  of 
A1  I. 


Table  5.  Heads  appearings  in  the  arc  spectrum  in  pure  Nt. 


kohi  ^ 

Int. 

i>ob$  cm-* 

reference 

i 

8741 

8 

11437 

violet  shade,  overlaid  by  Mn  i 
8740.93 

8522 

3 

11731 

violet  shade 

7821 

8 

12783 

red  shade 

7806 

3 

12808 

red  shade 

7743 

5 

12911 

violet  shade 

!  7347 

3 

13607 

red  shade 

7299 

3 

13696 

red  shade 

7150 

1 

13982 

violet  shade 

7040 

1 

14201 

red  shade 
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not  in  pure  nitrogen.  However  they  are  left  unanalysed;  overlapping  of  many 
bands  made  the  analysis  difiicvilt. 
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Aaatrmct 


Approximate  relations  between  C.  1.  E.  tristimulus  values  {X,  Y  and  Z)  of 
color  and  colar  temperature  (T)  of  Planckian  radiation  are  investgated  with 
numerous  colors.  Following  relations  are  found. 

X= a, + X  10«)  +  Cx(l/T  X  10*)2 

y=a»+Vl/7'xlO*)+C»(l/rxlO«)2  when  7= 3000—7000° K 

a,  +  6,  (1/7  X  10«) + C.  (1/7  X  10«)2 

For  these  empirical  equations,  some  theoretical  considerations  are  made. 

These  quadratic  forms  can  be  used  in  the  calculation  of  tristimulus  values 
of  color  of  Planckian  radiation  at  any  color  temperature. 


1.  Introduction 


Color  stimulus  values  of  an  object  are  changed  when  color  temperature  of 
the  illuminating  light  is  changed.  For  example,  even  if  the  tristimulus  values  of 
two  objects  are  equal  under  irradition  at  a  certain  color  temperature,  they  may 


b-JlO  • 
c  -HO  ' 
d  JSO  ' 
-e:2V’ 
f'-iX" 


Values  of  x 


Values  of  x 


Fig.  1.  Chromaticity  diagram. 
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be  different  under  irradiation  at  other  color  temperatures*. 

Thus,  as  shown  in  Figs,  la  and  lb,  various  colors  have  various  loci  with  the 
change  of  color  temperature  of  the  illuminating  light,  which  is  proved  in  actual 
cases.  For  example,  although  an  object  is  colored  under  a  certain  standard  light 
to  have  the  same  color  stimulus  values  as  those  of  a  given  sample  and  is  color- 
matched,  if  the  color  temperature  of  the  light  is  altered,  the  colors  may  no 
longer  match. 

The  relation  of  tristimulus  values  of  color  to  the  color  temperature  of  illum¬ 
inating  light  was  studied  and  empirical  formulas  which  would  be  useful  for  practical 
calculations  are  presented  in  this  work. 

2.  Theroretical  consideration 

2. 1.  Tristimulus  values  of  light  source  color. 

Tristimulus  values  of  light  source  color  are  defined  by 

L=^’‘^ErSi-pM  (1) 

where  Ei  is  the  radiation  energy  of  the  black  body  given  by  Placck’s  formula 
and  can  be  expanded  as  follows. 

. T 

=c./c..r/x{i-|;^+4(-§-)‘-  4(.§-)'+ . ]  (2 ) 

Si  is  an  approximate  function  for  tristimulus  values  of  equi-enengy  spectrum 
that  are  distribution  coeflScients  (x,  y,  z).  Peason’s  distribution  function  is  used 
for  Si,  which  is  defined  by 

S,=Ali^-exp[-qli]  {3) 

where  i  is  the  wave-length.  A,  P  and  q  are  constants. 

pi  in  (1)  is  the  spectral  characteristic  of  color  and  can  be  expressed  by  a  power 
series  expansion  of  spectral  reflection  or  transmission  curve  as 

pi  =  'SiBmA’^  (4) 

m-O 

therefore,  from  (1),  (2),  (3)  and  (4),  we  obtain 
L= J*  £1  •  S]  •  ^irfx = i4Ci(-^)  - 1(-^)  J 

*  In  this  discussion,  the  light  source  is  assumed  to  have  the  spectral  energy  distribution 
approximtely  equal  to  that  of  the  black  body. 
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putting  l!X=x 


-  r(F-m+3)  1/ 

C2\ 

V  c*  y 

L  2  \ 

T  ) 

nP-m+A)  ,  VCeV  r(P-m+5)  _  \  (  C2  \*  /’(P-w+6) 
12V  r  y 


720 


(^) 


(F— m+3)(P— w+4)  1 


(F— w+3X^— »»+4)(P— w+5)(F— w  +  6) 


720  \  T  / 

q* 

1  ( 

'  Ct  V  (F-m+3) 

•.••(F-mH-8) 

30240' 

.  r  ) 

9* 

1209600  v  r ) 


C2  V  (F-ni+3)----(F-w  +  10) 


(6) 


The  constants  F  and  9  in  (3)  are  given  by  P.  Moon  and  D. E.  Spencer'*-*’ 
as  shown  in  Table  1.  As  seen  in  this  table,  any  set  of  F  and  q  can  be  used  as 
the  approximate  values  of  spectral  distribution  coefficients  within  a  particular 
spectral  region.  Accordingly  the  whole  spectral  region  is  divided  into  several 
intervals. 

Table  1  Values  of  the  constants  in  Eq  (3) 


x^i)=xa  <i)  —XBit)-\-xcii):xA  (2)  570 

y{i)  =  yAH),  X^O.620  ft: y Ail)  — yau),  0.620^X^0.670  fi,  yn)  =  ycii),  X ^0.670  ft 


2<2) 

is  (i)  1  iCil) 

i/a  (2) 

VBH) 

Sc  (2) 

P 

Q 

365.3388 

164.9506 

500.0000  336.0385 

237.5000  199.1054 

182.1905 

100.9370 

3.2050 

0.12820 

1013.50 

679.045 

Let  us  use  just  one  appropriate  set  of  F  and  q  in  the  following  discussion. 
Now,  assuming  F— m>180and  T =3000~7000°K  (6)  can  be  written  approximately 
as  follows. 


1)  J.  Opt.  Soc.  Amer.,  33,  89-103  (1943). 

2)  J.  Opt.  Soc.  Amer.,  35,  399-427  (1945). 
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1  (  C^P-m+Z)  W  \  (  C*(/*-m+3)  V 
720\  Tq  /■^30240\  Tq  ) 


_ Lf. 

1209600  V 


Ci(P-m+3)  V 


Tq 


)-•] 


L=f,A-B^  Cx 


1 


C^P—m-\-Z) 

Tq 


]-' 


(7) 


Since  C*(/*— w+3)/7'9>4~5  from  Table  1 

L=±A-B^-Cx  -exi^-  ]  (8) 

As  the  result,  we  obtain  the  color  stimulus  value  Lo  for  a  bare  light  source 


L.=\'eSM=A-C,-  '~<P+3)  ■exp[--<^-<^+3)  j  ,3, 

Thus,  when  a  light  source  is  not  influenced  by  the  object  color,  l/T  dependence 
of  the  color  stimulus  value  is  approximately  represented  by  an  exponential  function 
at  least  under  the  assumed  conditions. 

2. 2  Tristimulus  values  of  object  color. 

Tristimulus  values  X,  Y,  Z,  of  an  object  are  defined  as  the  ratios  of  L  to 
the  stimulus  value  of  the  light  source  To. 


\  EiSipidX  \  Ei^iOioi  (  EiZxpidk 

- :  Y=Yn - ;  -  (10) 

^^EiVidX  ]^EiyidX  ^^ExyidX 

Since  coefficients  of  y  are  equal  in  (8)  and  (9),  Y  of  (10)  can  be  written  as  follows. 

Y =Bf,+Bx  exp  [Ci/T^vl+fix  exp  {2C2lTqy]+Bi  exp  [ZC^|Tqy\^ -  (11) 

which  can  be  expanded  as  below,  for  expiCzJTqy)  in  (11)  is  nearly  equal  to  unity 
and  m  is  small.  Thus  by  taking  the  first  terms  only,  we  have 

Y  ^(Bo-i-Bx+Bt+Bt+  •  •  •  Oh — ^^D{Bx+2B2+ZBt+  •  •  •  •) 

+  -^-^)V(Bi+4B,+9B,+  ----)  (12) 

where  C2iqy=D. 

From  these  considerations,  it  is  expected  that,  the  tristimulus  values  of  reflected 
color  of  an  object  of  known  spectral  reflection  characteristics  show  ayproximately 
linear  dependence  upon  reciprocal  color  temperature. 


Table  2.  Values  of  the  constants  in  empirical  equations 
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P.  B.  Blue  3^=  Polar  Brilliant  Blue  iR  160:  E.  Green  B  Supra  =  Erio  Green  B  Supra 
O.  Violet  3^=C)rasol  Violet  3/?:  O.  Blut  2GlV=Orasol  Blue  2^1^ 

O.  Orange  G/?lP=Orasol  Orange  GRW:  F.  Orange  G=Fat  Orange  G 


Table  3.  Chromatioity  coodinates  for  the  colors  in  standard  illuminant  (C) 
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P.  B.  Blue  3/?=  Polar  Briliant  Blue  3R  160%:  E.  Green  B  Supra  =  Erio  Green  B  Supra 
O.  Violet  3f?=C)rasol  Violet  3^:  O.  Orange  GRW  =Orasol  Orange  GRU' 


Percent  Reflectance 


Wavelength  in 
(g) 

Fig.  2.  Spectral  transmission  and  reflection  curves  of  samples. 
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For  Z  stimulus  value. 


qP  m^s  q^^* 


Hence  (10)  can  be  written  as 

Z = exp  [C. { q,{Pt + 3) - q^P, +3)) IT q^,[  (ko+ki  exp  [C2IT q,]  +  k2  exp  [2C2/r <7,] 

+  *»exp[3C2/7'«7.]+ -  Cl3) 

Since  exp(C*{<7,(P,+3)-^»(P.+3)}/T9»fl.l>exp[C2/T9»]>exp  [CtiTq,] 
the  factor  exp[Ci{q,iPp+3)—q^P,+3))ITq/},]  is  considerably  larger  than  the  factor 
exp[CtlTq,]  which  is  closer  to  unity  than  exp[C2/T?»]  in  (11).  Equations  (11)  and 
(13)  are  similar  in  forms  but  with  different  coefficient  for  different  object  colors. 

Hence,  the  general  feature  of  the  curves  1/T  vs  Z  stimulus  values  given  by 
equation  (13)  are  not  appreciably  altered  by  object  colors. 

As  for  X  stimulus  values,  distribution  curve  of  x  is  not  as  simple  as  that  of 
y  or  z,  but  as  for  their  relation  to  1/T,  it  is  similar  to  the  case  of  Z  stimulus 
values. 


3.  Empirical  equation  deduced  from  theoretical  consideration 

In  practice,  it  is  rather  difficult,  as  shown  in  the  preceding  section,  to  use  the 
theoretical  equations  for  calculation.  Therefore  it  will  be  desirable  to  formulate 
empirically  the  relation  between  object  color  and  reciprocal  color  temperature 
(mired)*  in  simple  rejjresentations.  Preceding  theoretical  condsideration  yields  such 
empirical  equations  which  will  be  shown  in  the  following  sections. 

The  results  obtained  by  numerical  calculations  are  shown  in  Table  2.  30 
selected  coordinate  method  was  used  which  is  most  usual  in  the  calculation  of 
tristimulus  values  of  color.  “  Standards  of  Color  ”  (made  by  Japxm  Color  Research 
Institute)  were  used  for  obtaining  the  spectral  reflection  curves,  and  dye-colored 
filters  for  the  spectral  transmission  curves. 

As  the  sample  numbers  of  color  charts  in  Table  2,  those  of  J.  C.  R.  I.  are  adopted. 

3. 1.  The  relation  between  Y  stimulus  value  (luminosity)  of  color  and  mired. 

As  seen  from  (12),  the  relation  between  Y  stimulus  value  and  mired  is  linear 
if  m  is  small,  that  is,  if  the  spectral  reflection  and  transmission  curves  of  this 
color  are  rather  simple. 

For  neutral  filter,  the  spectral  transmittance  and  reflectance  are  equal  to  Y 
stimulus  value.  Hence,  Y=Bo. 

For  various  samples,  numerical  calculations  were  performed  to  obtain  empirical 
equations,  which  are  shown  in  Table  2.  The  chromatid ty  coordinate  and  lum- 


*  Micro  reciprol  degree 
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inosity  of  the  samples  under  irradiation  of  standard  illuminant  “  C  ”  are  shown 
in  Table  3. 

In  Table  2,  standard  deviations  are  given  in  the  last  column.  Fig.  3  shows 
the  relation  between  tristimulus  values  and  mired. 

Linear  dependence  of  luminosity  upon  mired  was  verified  for  several  special 
samples. 

For  reflected  color  such  as  that  given  in  color  charts,  this  linearity  was  as¬ 
certained  for  all  the  samples. 

Dotted  curve  in  Fig.  2g  is  the  one  that  was  calculated  by  assuming  m  to  be 
large.  This  is  represented  by  a  quadratic  equation. 

3. 2.  The  relation  between  X,  Z  stimulus  values  of  color  and  mired. 

We  have  noted  before  that  the  relation  between  Z  stimulus  value  and  mired, 
(13),  is  not  appreciably  altered  by  the  difference  colors  of  objects.  Hence,  in  a 
limited  range  of  color  temperature,  the  quadratic  forms  can  be  used  with  a  suf¬ 
ficient  degree  of  accuracy.  They  are  shown  in  Table  2. 

Although  X  distribution  coefficients  in  the  X,  F,  Z  system  are  more  complicated 
in  this  case,  the  relations  between  X  stimulus  value  and  mired  are  similar  to  the 
case  of  Z  stimulus  value. 

For  example,  for  a  blue  or  green  colored  object,  x  distribution  coefficients  are 
predominantly  effective  only  in  a  short  wave  length  region,  and  for  red  or  orange 
object,  only  in  a  long  wave  length  region. 

For  the  intermediate  color,  though  x  distribution  coefficients  must  be  taken 
into  account  for  the  whole  visible  wave  length  region,  substantial  alteration  should 
not  occur  in  the  relation  between  X  stimulus  value  and  mired. 

Therefore  the  quadratic  formulae  are  also  applicable,  which  are  also  shown 
in  Table  2. 

3.3.  Empirical  equations  and  their  applications. 

As  shown  in  Table  2,  the  relation  between  mired  and  tristimulus  values  of 
various  colors  are  represented  by  the  following  empirical  equations: 


X =ax-^bxM 
Y =ay+byM  ->rCyM'^ 
Z =at+b,M +CtM^ 


Ox-h  bxM  -t-  CxM^ 

Oy  +  byM  -|- CyM* 
^  A+BM+CM^ 


where  coefficients  a,  b  and  c  are  characteristics  of  each  of  the  colors,  since 
A/=1/Txl0«,  A=ax+ay+a.,  B=bx+by+b,,  C=Cx+Cy+c,.  Applying  these  equations, 
color  stimulus  values  and  chromaticity  coordinates  are  easily  calculated  for  any 
color  temperature. 

In  practical  calculations,  tristimulus  values  for  “  A  ”,  “  B  ”  and  “  C  ”  standard 


On  the  Relation  Between  Tristimulus  Values  of  Color 


123 


illuminant  sources  are  numerically  calculated  at  first  by  the  usual  30  selected 
coordinate  method  using  colorimetric  computer,  thereby  the  coefficients  are 
determined. 

4.  Summary 

The  relation  of  tristimulus  values  of  color  and  color  temperature  of  a  light 
source  was  discussed  with  the  conclusion,  that  (in  the  color  temperature  region 
between  "A”  and  B"  standard  illuminant  sources)  this  relation  is  represented 
by  an  exponential  form  which  was  verified  for  various  color  samples  at  least 
within  the  accuracy  of  numerical  calculation. 

By  determining  tristimulus  values  of  an  object  color  for  "A”,  “B”  and 
“C”  standard  illuminant  sources  by  the  empirical  equation  obtained  and  then  by 
calculating  empirically  coefficients  for  this  color,  color  stimulus  values  can  be 
obtained  for  any  color. 
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